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Abstract  We  review  some  techniques  employed  in  the  studies  of  torques  due  to  bodily 
tides,  and  explain  why  the  MacDonald  formula  for  the  tidal  torque  is  valid  only  in  the  zeroth 
order  of  the  eccentricity  divided  by  the  quality  factor,  while  its  time-average  is  valid  in  the 
first  order.  As  a  result,  the  formula  cannot  be  used  for  analysis  in  higher  orders  of  e/ Q. 
This  necessitates  some  corrections  in  the  current  theory  of  tidal  despinning  and  libration 
damping  (though  the  qualitative  conclusions  of  that  theory  may  largely  remain  correct).  We 
demonstrate  that  in  the  case  when  the  inclinations  are  small  and  the  phase  lags  of  the  tidal 
harmonics  are  proportional  to  the  frequency,  the  Darwin-Kaula  expansion  is  equivalent  to 
a  corrected  version  of  the  MacDonald  method.  The  latter  method  rests  on  the  assumption 
of  existence  of  one  total  double  bulge.  The  necessary  correction  to  MacDonald’s  approach 
would  be  to  assert  (following  Singer,  Geophys.  J.  R.  Astron.  Soc.,  15:  205-226,  1968)  that 
the  phase  lag  of  this  integral  bulge  is  not  constant,  but  is  proportional  to  the  instantaneous 
synodal  frequency  (which  is  twice  the  difference  between  the  evolution  rates  of  the  true 
anomaly  and  the  sidereal  angle).  This  equivalence  of  two  descriptions  becomes  violated  by  a 
nonlinear  dependence  of  the  phase  lag  upon  the  tidal  frequency.  It  remains  unclear  whether  it 
is  violated  at  higher  inclinations.  Another  goal  of  our  paper  is  to  compare  two  derivations  of  a 
popular  formula  for  the  tidal  despinning  rate,  and  emphasise  that  both  are  strongly  limited  to 
the  case  of  a  vanishing  inclination  and  a  certain  (sadly,  unrealistic)  law  of  frequency-depen¬ 
dence  of  the  quality  factor  Q — the  law  that  follows  from  the  phase  lag  being  proportional  to 
frequency.  One  of  the  said  derivations  is  based  on  the  MacDonald  torque,  the  other  on  the 
Darwin  torque.  Fortunately,  the  second  approach  is  general  enough  to  accommodate  both  a 
finite  inclination  and  the  actual  rheology.  We  also  address  the  rheological  models  with  the  Q 
factor  scaling  as  the  tidal  frequency  to  a  positive  fractional  power,  and  disprove  the  popular 
belief  that  these  models  introduce  discontinuities  into  the  equations  and  thus  are  unrealistic 
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at  low  frequencies.  Although  such  models  indeed  make  the  conventional  expressions  for  the 
torque  diverge  at  vanishing  frequencies,  the  emerging  infinities  reveal  not  the  impossible 
nature  of  one  or  another  rheology,  but  a  subtle  flaw  in  the  underlying  mathematical  model 
of  friction.  Flawed  is  the  common  misassumption  that  damping  merely  provides  phase  lags 
to  the  terms  of  the  Fourier  series  for  the  tidal  potential.  A  careful  hydrodynamical  treatment 
by  Sir  George  Darwin  (1879),  with  viscosity  explicitly  included,  had  demonstrated  that  the 
magnitudes  of  the  terms,  too,  get  changed — a  fine  detail  later  neglected  as  “irrelevant”.  Rein¬ 
stating  of  this  detail  tames  the  fake  infinities  and  rehabilitates  the  “impossible”  scaling  law 
(which  happens  to  be  the  actual  law  the  terrestrial  planets  obey  at  low  frequencies).  Finally, 
we  explore  the  limitations  of  the  popular  formula  interconnecting  the  quality  factor  and  the 
phase  lag.  It  turns  out  that,  for  low  values  of  Q,  the  quality  factor  is  no  longer  equal  to  the 
cotangent  of  the  lag. 

Keywords  Tides  •  Body  tides  •  Bodily  tides  ■  Land  tides  •  Tidal  forces  •  Tidal  torques  • 
MacDonald  torques  ■  Libration  •  Natural  satellites  •  Tidal  despinning  •  Spin-orbit  interaction  • 
Spin-orbit  coupling  •  Spin-orbit  resonances 


1  Prologue 


When  it  shall  be  found  that  much  is  omitted, 
let  it  not  be  forgotten  that  much  likewise  is  performed 

Samuel  Johnson,  1755 

In  his  short  work  “Untersuchung  der  Frage  known  among  the  historians  also  as  the 
“ Spin-Cycle  essay,”  Immanuel  Kant  (1754)  stated  that  the  Moon  not  only  pulls  the  Earth,  but 
also  exerts  a  retarding  torque  upon  its  surface;  this  torque  slows  down  the  Earth’s  rotation 
and  lets  go  only  when  terrestrial  days  become  as  long  as  lunar  months.  Although  Kant  had  in 
mind  only  the  ocean  tides,  not  the  bodily  ones,  we  may  say  that,  qualitatively,  he  predicted  the 
celebrated  1 : 1  spin-orbit  resonance,  the  pas  de  deux  wherein  Pluto  and  Charon  are  locked. 

For  the  first  time,  the  idea  of  tidal  action  not  being  confined  only  to  the  fluid  portion  of  the 
planet  but  affecting  also  the  solid,  so  as  to  induce  a  state  of  varying  strain,  was  put  forward 
by  John  Herschel  (son  of  astronomer  William  Herschel),  as  a  minor  aside  in  a  paper  devoted 
to  volcanism  and  earthquakes  (Herschel  1863).  The  earliest  mathematical  description  of  land 
tides  in  their  dynamics  was  offered  by  George  Darwin  (son  of  naturalist  Charles  Darwin  and 
great-grandson  of  poet  and  philosopher  Erasmus  Darwin). 

Following  his  predecessors  Roche  (1849)  and  Thomson  (1863),  who  had  calculated  the 
figure  of  a  static  tide,  Darwin  (1879)  assumed  the  Earth  homogeneous  and  consisting  of 
an  incompressible  fluid.  To  account  for  dynamics,  he  also  assumed  that  the  viscosity  was 
the  sole  source  of  the  tidal  friction.  Relying  on  this  model,  Darwin  (1880,  1908)  derived  a 
tide-generated  disturbing  potential  expanded  into  a  Fourier  series.  Substitution  thereof  into 
the  Lagrange-type  planetary  equations  led  him  to  expressions  for  the  time  derivatives  of  the 
orbital  elements  via  partial  derivatives  of  the  disturbing  potential  with  respect  to  the  elements. 

An  impressive  generalisation  of  Darwin’s  work  by  Kaula  (1964),  and  the  subsequent 
flow  of  new  concepts  and  applications  (MacDonald  1964;  Goldreich  1966a,b;  Goldreich  and 
Peale  1996;  Singer  1968;  Mignard  1979,  1980;  Touma  and  Wisdom  1994;  Neron  de  Sur- 
gy  and  Laskar  1997;  Krasinsky  2002,  2006;  Getino  et  al.  2003;  Ferraz-Mello  et  al.  2008; 
Efroimsky  2008)  made  bodily  tides  a  rapidly  developing  area  of  the  planetary  astronomy. 
The  vast  and  growing  volume  of  the  relevant  material  leaves  one  no  chance  to  glean  it  all  in 
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one  review.  Therefore,  we  shall  concentrate  on  one  special  aspect  of  this  research,  the  tidal 
torques  emerging  from  the  bodily  tides.  Moreover,  we  shall  dwell  solely  on  the  techniques, 
not  applications. 

Although  our  review  will  at  times  be  very  critical,  it  should  from  the  beginning  be  agreed 
that  our  criticisms  are  intended  in  the  spirit  of  the  above  quotation  from  Samuel  Johnson. 

Along  with  reviewing  the  current  state  of  the  field,  we  shall  provide  some  new  results  of 
our  own.  Specifically,  we  shall  address  the  rheological  models  with  the  Q  factor  scaling  as 
the  tidal  frequency  to  a  positive  fractional  exponential.  We  shall  demonstrate  that,  contrary  to 
the  common  opinion,  such  rheologies  do  not  cause  infinities  in  the  expression  for  the  torque. 
We  shall  also  derive  an  expression  for  the  tidal  torque  decelerating  a  terrestrial  planet  obeying 
such  a  rheology.  (That  the  realistic  terrestrial  bodies  indeed  obey  this  class  of  rheologies  has 
been  explained  in  Efroimsky  and  Lainey  (2007).) 


2  Trivia 


In  this  section,  we  shall  briefly  recall  how  a  satellite-generated  potential  in  a  point  on  or 
inside  the  planet  is  expressed  through  the  latitude,  longitude,  and  the  radial  distance  of  the 
point. 

Let  us  begin  from  the  first  principles.  The  dynamics  of  point  masses  m,  located  at  iner¬ 
tial-frame-related  positions  jo, , 


mtPj  =  mi  y  .  Gm.j  — ,  rtj  =  pj  -  pt,  i,  j  =  1, . . .,  N, 
j¥=‘  r0 

may  be  conveniently  reformulated  in  terms  of  the  relative-to-the-primary  locations 

ri  =  r0i  =  p,  -  p0, 

p0  standing  for  the  position  of  the  primary.  The  difference  between 

mo?  ,o 


j^i.O  rij 


'i  0 


and 


amounts  to: 


Jr0j  .  m/ /■(),' 


*0=2 


j#,  o  o  J 


r3 
r0  i 


h=  X  G^-  X  g'^-G 
J¥=‘.  o 

Uj  being  the  potential: 


mjfj  (ttii  +  mo)  r _  3(7; 


./' '  i.o  riJ  j  >  ;.o  ri 


d?i 


G  (mi  +  mo) 

Ui  = -  +  Wt, 


with  the  disturbance 


Z1  ?i  ■  r  i 

Gm  1 

i¥=i 


(1) 

(2) 

(3) 


(4) 


(5) 


(6) 

(7) 
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singled  out.  This  disturbing  potential  acting  on  mass  m,-  is  generated  by  the  masses  m  j  other 
than  m,  or  the  primary.  It  deviates  from  the  Newtonian  one  by  the  amendment  GmjrJJrj  ■  r  j 
emerging  in  the  noninertial  frame  associated  with  the  primary. 

In  the  simplest  case  of  one  secondary,  a  satellite  of  mass  mi  =  M*aV  located  at  a  plane- 
tocentric  position  r\  =  r* ,  will  be  creating  at  some  point  ?2  =  R  a  perturbing  potential 


W(R,  r)  =  -GM* 


R  r 


I R  -  r 


(8) 


expandable  over  the  Legendre  polynomials  (for  R  <  r*)  by  means  of  the  formulae 


1 

\R-F\ 


1  ^  /  R\' 

^rZ(^)  p'(cosy) 
1=0  v  7 


(9) 


and 


R  r 


Rr*  cos  y 


1  R 


-Pi  (cosy), 


(10) 


y  being  the  angular  separation  between  R  and  r* ,  subtended  at  the  point  of  origin,  which  we 
shall  naturally  choose  to  coincide  with  the  planet’s  center  of  mass.  Together,  the  former  and 
the  latter  formulae  yield: 


W(R,r)  =  -- 


GM;at  y(R\ 

r*  1  \  r*  I 
1=2  v  7 


—  )  P,(cosy), 


(11) 


where  we  have  neglected  the  1  =  0  term  —GM*at/r* ,  because  it  bears  no  dependence  upon 
R,  and  in  practical  problems  is  attributed  to  the  principal  part  of  the  potential,  not  to  the  one 
regarded  as  perturbation.  The  angle  y  can  be  expressed  via  spherical  coordinates  as: 


cos  y  = 


R  r 


Rr* 


sin  (f>  sin  <p*  +  cos  cf>  cos  <f>*  cos(L  —  L*), 


(12) 


(R,  <t>,  X)  being  the  planetocentric  distance,  the  latitude,  and  the  longitude  of  the  point  where 
the  disturbance  is  experienced;  and  (r* ,  <p* ,  k*)  being  the  spherical  coordinates  of  the  satellite. 
It  is  customary  (though  not  at  all  obligatory)  to  reckon  the  longitudes  from  a  planet-fixed 
meridian,  in  which  case  the  subsequent  formulae  for  the  potential  come  out  written  in  a 
reference  frame  co-rotating  with  the  planet. 

A  Legendre  polynomial  of  cos  y,  too,  can  be  expressed  via  the  spherical  coordinates: 


'  n  _  m\\ 

Pi  (cosy)  =  ^  n  |  ^’(2  -  80m)  Pi  m  (sin  cf>)Pim  (sin  cp*)  cos  m(k  -  k*), 


m=  0 


(/  +  m)\ 


(13) 


substitution  whereof  into  (11)  results  in 

GM: 


W(R,  r  )  =  -- 


1=2  V  7  m=  0 


(Z  —  m)\ 
(l  +  m)\ 


(2  —  <5o,„) 


x  Pi m (sin  (p)P]m (sin  <p*)  cos  m(k  —  k*). 


(14) 


Evidently,  this  formalism  will  stay  unaltered,  if  the  role  of  the  tide-raising  satellite  is  played 
by  the  Sun,  or  by  another  satellite,  or  by  another  planet.  (In  this  case,  what  we  call  Msat  will,  in 
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fact,  denote  the  mass  of  the  Sun,  or  of  the  other  satellite,  or  of  the  other  planet.)  Likewise,  the 
formalism  may  in  its  entirety  be  applied  to  a  satellite  regarded  as  a  tidally-disturbed  primary, 
the  planet  being  treated  as  a  tide-raising  body  (and  M*al  now  standing  for  the  planetary  mass). 


3  The  Kaula  expansion  for  a  tidal  potential 


Kaula  (1961)  came  up  with  a  remarkable  formula 

H-i  _  1  /  i  \?+i  00 

P/(sin</>*)  cosmk*  +  sf—  1  sinmk*  =  I  —  |  ^  TjmpO*) 

' a  '  p= o 


OO 

Z  G^e*) 

q= — oo 


cos  (<mpq  -  m6*)  +  sin  (<mpq  -  me*) 

Sin  (<mpq  “  m6*)  ~  C°S  (<mpq  “  m9*) 


l-m  even 


(15) 


l—m  odd 


where  Fimp(i)  are  the  inclination  functions  (Gooding  2008);  Gipq(e)  are  the  eccentricity 
polynomials  identical  to  the  Hansen  coefficients  ;  the  notation  T  is  used  to 

avoid  confusion  with  the  inclination;  and  the  auxiliary  combinations  i>*mpq  are  defined  as:1 

Hmpq  =  d-  2 p)a>*  +  (l  -  2 p  +  q)M*  +  mQ*.  (16) 

This  development  enabled  Kaula  (1961,  1964)  to  carry  out  a  transformation  from  the  tide¬ 
raising  satellite’s  spherical  coordinates  to  its  orbital  elements  and  the  sidereal  time  6*.  These 
elements  (the  semimajor  axis  a*,  eccentricity  e* ,  inclination  i* ,  periapse  co* ,  ascending  node 
£2*,  mean  anomaly  A4*)  are  introduced  in  a  frame  that  is  associated  with  the  equator  but  is 
not  co-rotating  with  it.  In  terms  of  these  parameters. 


W(R,r*) 


gm*m 

a*  ^\a* ) 

1=2  v  7  m= 0 


(l  —  m)\ 
(, l  +  m)\ 


(2  -  Som)Ph 


,(sin</>)  z  Plmp(1*) 

P= 0 


oo  r  f  l  l~m  even 

Z  GIpq(e*)  cosml  “S  Kmpq  ~ md*) 

q~^oo  L  I  sln  I  l—m  odd  V  ’ 


l—m  even 

+  sin  m  X 

sin 

—  cos 

(<npq  -  ™e*) 

\—m  odd 

or,  after  carrying  out  the  multiplication  of  the  sine  and  cosine  functions: 


W(R,r)  =  -- 


G  -y  /  R  \  yo 
a*  \a*J 

7=2  v  7  m= 0 


(/  —  m)\ 
(/  +  m)\ 


(2  -  t$0m)  Plm  (sin  <p) 


Z^lmpO'*)  Z  Glpq(fi*) 

p= 0  q=—o o 


COS 

sin 


l—m  even 
l— m  odd 


(wSnpq-m^  +  O)  (18) 


1  This  definition  agrees  with  that  by  Kaula  (1961,  1964,  1966),  but  differs  from  the  one  by  Lambeck  (1980) 
who  incorporated  —m6*  into  uj^pq- 
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4  Physical  assumptions  involved  in  Kaula’s  theory 

If  the  primary  is  not  a  point  mass,  it  becomes  distorted  by  potential  W(R,  ?*).  The  distor¬ 
tion  of  shape  will,  in  its  turn,  generate  some  extra  potential  perturbation  whose  calculation  is 
complicated  by  the  tide-raising  potential  (18)  evolving  in  time  and  having  a  rich  spectrum  of 
frequencies.  The  response  of  the  primary’s  shape  to  each  of  these  is  different  and  depends  on 
the  properties  of  the  planet’s  material.  This  is  a  situation  where  the  linear  approach  becomes 
most  helpful,  when  applicable.2 

The  linear  theory  of  bodily  tides  comprises  two  independent  assertions.  One  is  that  the 
energy  attenuation  rate  (E)  at  each  harmonic  depends  solely  on  the  frequency  x  and  on  the 
amplitude  Speak  (x)>  and  is  not  influenced  by  the  rest  of  the  spectrum.  This  is  written  down  as 
(£(x))  =  -x£Peak(x)/£Kx)>  which  is  equivalent  to  A£cycie(x)  =  -2jr£peak(x)/G(x)> 
where  A  Ecyc\e  (x )  is  the  one-cycle  energy  loss,  and  Q  (x )  is  the  quality  factor.  The  other  asser¬ 
tion  is  that  each  stationary  tidal  change  of  the  potential,  W/,  inflicts  on  the  planet’s  shape  a  lin¬ 
ear  deformation.  Each  of  these  deformations,  in  their  turn,  amend  the  potential  of  the  primary 
with  an  addition  proportional  to  the  Love  number  k /.  As  known  from  the  potential  theory,  an 
addition  proportional  to  Pj  (cos  y )  must  be  decreasing  outside  the  spherical  primary  as  r  ~ (  +  ^ . 
Hence,  were  the  external  potential  perturbation  W  static  (or,  equivalently,  were  the  response 
of  the  material  instant),  the  tidal  addition  to  the  planetary  potential  would  have  assumed  the 
form3 


U(r)  = 


W,(R,  r) 


™  /  n\  1+ 

5‘(‘) 

-£*(?) 

1=2  v  7 

I  yv  r  i  i—m  even 

2>tapG*>  X  ^P^*)  dn  Knpq -«(*  +  **))  (19) 

p= 0  q=~oc  I  Sln  I  1  -m  odd  V  2 


(D'i 

\  /  I'M — n 


(/  —  m)\ 
(/  +  m)\ 


=0 
/— m  even 


R  being  the  mean  equatorial  (equivolumetric)  radius  of  the  planet,  R  =  ( R ,  0,  X)  being  a 
particular  surface  point,  and  r  =  (r,  <f>,  X)  being  an  exterior  point  located  right  above  the 
surface  point  I?,  at  a  planetocentric  radius  r  >  R. 

As  we  intend  to  study  the  effect  of  this  potential  on  another  external  body,  a  similar  trans¬ 
formation  should  be  applied  to  the  coordinates  (r,  0,  X),  to  express  W  through  the  orbital 
elements  of  this  body.  Employment  of  (15),  this  time  not  for  r  but  for  r,  leads  to: 


2  For  most  materials,  departure  from  linearity  becomes  considerable  when  the  strains  approach  10  . 

(Karato  2007) 

3  Following  MacDonald  (1964)  and  Singer  (1968),  we  denote  the  tide-raising  potential  with  W  and  the 
bodily-tide  one  with  U .  In  his  original  paper,  Kaula  (1964)  called  these  potentials  U  and  T .  while  in  the 
book  he  switched  to  U  and  UT  (Kaula  1968).  Be  mindful  that  we  are  using  a  sign  convention  different 
from  that  of  Kaula.  As  our  forces  are  negative  gradients  of  potentials,  our  potentials  are  negative  to  those  of 
Kaula. 
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U(r)  =  - 


1=2  x  x  '  m= 0 


( l  —  m)\ 
(, l  +  m)\ 


(2  -  S0m) 


"y  7*lmp(i  )  ^  Gipq(e  )  y  '  FlmhlO 


P= 0 


q=—oo 


h=0 


oo 

X  Gihj(e)  cos  -  m6>*)  -  (i>imhj  -  m0)l , 


(20) 


j=-oo 


a  formula  that  generalises  the  tidal  theory  of  Darwin  (1908,  p.  334)  to  /  and  |^|  larger  than 
2.  Both  Kaula  (1964),  who  derived  this  milestone  result,  and  Darwin,  who  had  developed  its 
simplified  version,  realised  that  this  machinery  would  work  only  after  the  material’s  delayed 
reaction  to  perturbation  (18)  is  somehow  taken  into  account.  Until  then  (20)  remains  idealised, 
in  that  it  corresponds  to  an  unphysical  case  of  instantaneous  response. 

To  account  for  damping,  Kaula  (1964)  followed  the  path  of  Darwin  (1880,  1908):  he 
endowed  each  term  of  the  Fourier  series  (20)  with  a  real  phase  lag  of  its  own,  eimpq,  where¬ 
after  the  ultimate  form  of  Kaula’s  expansion  became 


U(r)  =  - 


1=2  X  7  X  7  m= 0 


(/  —  m)\ 


(2  -  (5om) 


(/  +  m)\ 

X  Glpq(e*)XFlmh(0 

p= 0  q=-oo  h=0 

oo 

y  \  Gihj(^)  COS  ^Uimpq  mG  ^  —  (r’lmhj  fH&)  ^lmpqj- 


(21) 


J=-o o 


This  empirical  method  of  including  dissipation  into  the  picture  contains  in  itself  an  important 
omission,  of  which  Sir  George  Darwin  was  aware,  but  which  was  overlooked  by  his  succes¬ 
sors.  Briefly  speaking,  even  in  a  linear  system  a  dissipation  process  is  not  fully  accounted  for 
by  amending  phases  of  the  Fourier  components.  This  observation  happens  to  be  of  relevance 
in  the  theory  of  tidal  torques.  We  shall  return  to  this  point  in  Sect.  9. 


5  The  two  sidereal  angles 

Kaula’s  construction  contains  a  seemingly  redundant  fixture,  which  turns  out  to  be  an  impor¬ 
tant  and  useful  acquisition.  This  is  Kaula’s  introducing  two  sidereal  angles  instead  of  one. 
As  these  angles,  6  and  9*,  are  not  orbital  elements  of  the  tide-raising  and  tidally  disturbed 
moons,  but  are  parameters  characterizing  the  instantaneous  attitude  of  the  planet,  it  may  look 
strange  that  Kaula  (1964)  assumed  them  to  be  different  entities.  To  understand  his  point,  let 
us  trace  the  physical  origin  of  the  phase  lag.  The  material  of  the  primary  is  being  deformed  by 
a  tidal  stress  whose  spectrum  contains  an  infinite  number  of  frequencies,  the  reaction  of  the 
material  to  each  of  these  being  different.  In  a  linear  regime,  the  strain  has  the  same  spectrum, 
with  each  harmonic  delayed  by  its  own  time  lag  Aqmpq.  Singer  (1968),  and  later  Mignard 
(1979,  1980),  assumed  that  all  Arimpq  are  equal  to  one  another:  Afimpq  =  At.  If  this  were 
true,  then  in  Kaula’s  series  each  argument  —  mG*  would  have  to  be  substituted  with 
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^(delayed) 

^lmpq 


^(delayed) 

—  mO 


tlfmpqO  -  At)  -  mO*(t  -  At) 

«4pq(0  -  me*(t)  -  [ffmpq  -  md*]  At 

ftapqCO  -  mO*(t)  -  [(/  -  2 p)w*  +  (l  -  2p  +  q)M* 

+  m(Q.*  -  0*)]  At.  (22) 


In  reality,  however,  the  time  lag  is  a  function  of  frequency,  for  which  reason  the  delays 
Afimpq  will  be  different  for  each  harmonic  involved.  This  is  why  the  arguments  u1*lpq  —  mO* 
at  the  moment  t  should  rather  be  replaced  with 


^(delayed) 

^lmpq 


^  (delayed) 
~~  m^lmpq 


<Pq  -mG*-[(l-2p)d>* 

+  (/  —  2 p  +  +  TTl(£ 2*  —  $*)]  A/]mpq. 


Specifically, 

utapqyed)  =  ^tmpq  “  [(l  ~  2P )«*  +  U  ~  2 P  +  q)M*  +  mf2*]  Arimpq 

and 


(delayed)  _  • 

^lmpq  —  “  &  A'Impq! 

9*  being  the  planet  spin  rate.  In  brief,  (23)  can  be  rewritten  as 

^(delayed)  (delayed)  *  * 

l'lmpq  —  W^lmpq  =  ’Vpq  —  mt>  “  <ulmpqAtlmpq. 

We  see  that  the  total  phase  lags  eimpq  introduced  by  Kaula  are  given  by 

eimpq  =  [(/  -  2 p)w*  +  (Z  -  2p  +  q)M*  +  m{Q.*  -  0*)]  Aqmpq 

=  ^Impq  Ahmpq  =  ±Xlmpq  ^fimpq, 

the  tidal  harmonic  &>*mpq  being  introduced  as 

wfmpq  =  (/  -  2 P)a*  +  (l-2p  +  q)M*  +  m(Q*  -  0*), 
the  positively-defined  physical  frequency 

Xlmpq  =  Kpql  =  I (l-  2 p)a>*  +  (l -2p  +  q)M*  +  m{Q*  -  9*)\ 


(23) 


(24) 


(25) 

(26) 


being  the  actual  physical  Impq  tidal  frequency  excited  in  the  primary’s  material.  The  appro¬ 
priate  positively-defined  time  delay  Afimpq  depends  on  this  physical  frequency,  for  which 
reason  the  delays  Afimpq  are,  generally,  different  from  one  another.4 * 

The  sign  on  the  right-hand  side  of  (24)  is  simply  the  sign  of  &)*mpq.  The  sign  evidently 
depends  on  whether  mO  falls  short  of  or  exceeds  the  linear  combination  (Z  —  2 p)d>*  + 
(l-2p  +  q)M*  +  mQ.*  «  (Z  -  2p  +  q)M*. 

The  origin  and  meaning  of  the  phase  lag  6impq  being  now  transparent,  one  may  express 
the  cosine  functions  in  (21)  either  as 


'[(^lmpq  m@  )  (ulmhj  m@)  6lmpqJ 


(27) 


4  When  Kaula  was  developing  his  theory,  the  functional  form  of  the  dependence  A t(x)  was  not  yet  known. 

Reliable  data  became  available  only  in  the  final  quarter  of  the  past  century.  See  formula  (83)  below. 
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(where  9*  and  0  are  identical  and  cancel  one  another),  or  simply  as 

cos  [(C*  -  -»<")  -  (»-)  - '”«)]  ■  rn 

In  (28)  we  have  the  delayed  siderial  angle,  separated  from  the  actual  angle,  9,  by 

-6  Afimpq,  the  time  lag  Aqmpq  being  a  function  of  ximpq  =  Nmpql- 


6  The  Darwin-Kaula-Goldreich  expansion  for  the  tidal  torque 


Now  we  are  prepared  to  calculate  the  planet-perturbing  tidal  torque.  Since  in  what  follows  we 
shall  dwell  on  the  low-inclination  case,  it  will  be  sufficient  to  derive  the  torque’s  component 
orthogonal  to  the  planetary  equator: 


r  = 


—  Afsat 


dU(r) 

99 


(29) 


Msat  being  the  mass  of  the  tide-disturbed  satellite,  and  the  “minus”  sign  emerging  due  to  our 
choice  not  of  the  astronomical  but  of  the  physical  sign  convention.  Adoption  of  the  latter 
convention  implies  the  emergence  of  a  “minus”  sign  in  the  expression  for  the  potential  of  a 
point  mass:  —GM/r.  This  “minus”  sign  then  shows  up  on  the  right-hand  sides  of  (6-8)  and, 
later,  of  (19-21).  It  is  then  compensated  by  the  “minus”  sign  standing  in  (29). 

The  right  way  of  calculating  9U ( r)/99  is  to  take  the  derivative  of  (28)  with  respect  to  9, 
then  to  insert  (23)  into  the  result,  and  finally  to  get  rid  of  the  sidereal  angle  completely,  by 
imposing  the  constraint  9*  =  9.  This  will  yield:5 


1=2  V  7  X  7  m=  0 

I  oo  I 

X^impO'*)  X  Glpq(e*)  l>tah(0 

p= 0  q=-o o  h= 0 


(/  —  m)\ 
(/  +  m)\ 


2m 


oo 

^  ,  ^lhj(^)  sin  ^lmpq  —  ftlmhj  —  ^lmpqj* 

j=~o o 


(30) 


In  the  case  of  the  tide-raising  satellite  coinciding  with  the  tide-perturbed  one,  Msat  =  M*at, 
and  all  the  elements  become  identical  to  their  counterparts  with  an  asterisk.  For  a  primary 
body  not  in  a  tidal  lock  with  its  satellite,6  the  torque  (30)  can  be  split  into  two  parts.  The 


5  Formally,  one  can  as  well  differentiate  (27)  instead  of  (28),  first  ignoring  the  fact  that  6 *  and  6  are  identical 
and  then,  after  differentiation,  permitting  them  to  cancel  one  another.  Though  this  method  produces  the  same 
result  as  the  rigorous  calculation,  it  nonetheless  remains  a  formal  procedure  lacking  physics  in  it. 

6  This  caveat  is  relevant,  because  in  resonances  expressions  (33-35)  will  require  modifications.  For  example, 
the  sidereal  angle  of  a  satellite  tidally  locked  in  a  1  :  1  resonance  will  be:  6  =  £2  +  co  +  Ai  +  1 80°  +  a  +  O  (z‘2), 
letter  a  denoting  the  librating  angle,  which  is  subject  to  damping  and  therefore  is  normally  small  (less  than 
2”  for  the  Moon).  Inserting  the  said  formula  for  6  into  the  expression  (25)  for  the  tidal  harmonic,  we  obtain, 
in  neglect  of  —mcr. 

"lmpq  =  (l-2p-  m)u*  +  (l-2p  +  q-  m)M*. 

We  now  see  that,  since  0  is  a  function  of  the  other  angles,  different  sets  of  the  indices’s  values  will  correspond 
to  one  value  of  the  tidal  frequency.  We  shall  illustrate  this  by  considering  the  so-called  anomalistic  modes  ±  A4 
in  the  potential.  These  modes,  corresponding  to  the  physical  frequency  \Ai\,  are  given  by  ( Impq )  =  (201 ,  d=l) 
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first  part  is  constituted  by  those  terms  of  (30),  in  which  indices  ( p ,  q )  coincide  with  (h,  j), 
and  therefore  all  uimhj  cancel  with  uj^pq,  provided  the  tidally-perturbed  satellite  and  the  tide¬ 
raising  one  are  the  same  body.  This  component  of  the  torque  is,  therefore,  constant.  The  rest 
of  the  total  sum  (30)  will  be  denoted  with  f .  It  is  comprised  of  the  terms,  in  which  the  pairs 
( p ,  q)  differ  from  (h,  j).  Accordingly,  these  terms  contain  the  differences 

ffmpq  -  Hmhj  =  (/  -  2 p  +  q)M*  -  (/  -  2 h  +  j)M  +  m{Q.*  -  £2) 

+  /(&)*  —  a>)  —  2pa>*  +  2hu>.  (31) 

When  the  tidally-perturbed  and  tide-raising  moons  are  the  same  body,  this  becomes 

njnipq  -  I'lmhj  =  (2/i  -  2p  +  q  -  j)M  +  (2 h  -  2 p)a>,  (32) 

whence  we  see  that  the  oscillating  component  of  the  torque,  f,  consists  of  two  parts.  The 
part  with  h  —  p  =  0  and  q  —  j  ^  0  consists  solely  of  short-period  terms,  and  it  averages  out 
trivially. 

The  mixed-period  part  of  (30),  with  h  —  p  ^  0,  consists  of  both  a  short-period  contribution 
dependent  upon  the  mean  anomaly,  and  a  long-period  contribution  depending  upon  the  argu¬ 
ment  of  pericentre.  All  such  terms  contain  multipliers  like  F2mp  (i) F2mh  (0,  where  h  p,  and 
F220  =  3  +  0(i* 2),  F2w  =  3/2 sin i  +  0(/3),  F2\\  =  -3/2 sin /+  0(i3),  F22\  =  3/2 sin2  i, 
the  other  relevant  F2mn ’s  being  of  higher  order  than  O  (r ) .  So  the  only  long-period  terms  that 
we  have  to  consider  in  (30)  involve  products:  //>io(0 F2\ 1  (0>  T211  (i)F2io(i),  T220 (0 ^221  (i)> 
and  F221  (0*220  (0-  However,  these  products  are  of  order  0(i2).  Thus,  while  the  short-period 
terms  in  (30)  average  out  over  one  rotation  period  of  the  moon  about  the  planet,  the  long- 
period  terms  are  of  order  0{e2i2),  the  e 2  coming  from  the  G2pq  functions.  (Indeed,  when  h 
and  p  differ  by  1,  then  q  and  j  must  differ  by  2,  to  eliminate  the  mean  anomaly,  i.e.,  to  make 
the  term  long-period  and  not  short-period.)  So  both  the  short-  and  long-period  contributions 
may  be  neglected  in  our  approximation.7 


Footnote  6  continued 

and  also  by  ( Impq )  =  (220,  ±1).  Although  the  m  =  0  terms  enter  the  potential,  they  will  not  be  in  the  torque, 

as  can  be  observed  by  differentiating  (19)  with  respect  to  X,  or  by  differentiating  equation  (21)  with  respect  to 
—6,  or  simply  by  noticing  the  presence  of  the  factor  m  on  the  right-hand  side  of  (30).  Nonetheless,  we  see  that 

there  exists  a  pair  of  m  =  2  terms,  which  provides  an  anomalistic  input  into  the  torque.  This  way,  the  case  of 

libration  deserves  a  separate  consideration,  as  it  is  more  involved  than  that  of  tidal  despinning.  Specifically, 

in  the  case  of  libration  a  value  of  the  tidal  frequency  may  correspond  to  different  sets  of  the  indices’  values. 

7  Had  we  tried  to  expand  our  treatment  to  higher  inclinations,  our  neglect  of  the  short-period  terms  would 
remain  legitimate,  for  they  still  would  average  out  over  one  rotation  period  of  the  satellite  about  its  primary. 
As  for  the  long-period  terms,  it  would  be  tempting  to  say  that  these  average  out  over  the  apsidal-precession 
period.  The  latter  is  much  shorter  than  the  time  scale  of  the  planetary  spin  deceleration,  a  circumstance  that 
may  seem  a  safe  justification  for  the  neglect  of  the  long-period  terms  also  for  higher  inclinations.  However,  a 
word  of  warning  would  be  appropriate  here.  As  well  known  from  Kozai  (1959a),  who  took  into  account  the 
primary’s  nonsphericity,  the  pericentre  of  a  satellite  inclined  by  about  63°  or  117°  will  neither  advance  nor 
retard,  at  least  within  the  first-order  (in  J2 )  perturbation  theory.  (For  a  critical  review  of  Kozai ’s  theory  see 
Taff  (1985).)  Kozai’s  original  attempt  to  introduce  corrections  owing  to  73  and  74  was  flawed  because  in  the 
vicinity  of  the  critical  inclinations  these  terms  should  be  considered  not  as  higher-order  but  rather  as  leading. 
His  later  analysis  demonstrated  that  at  these  inclinations  the  satellite’s  perigee  should  librate  about  90°  or 
270°  (Kozai  1962).  Under  these  circumstances,  the  long-period  terms  in  our  expression  for  the  torque  will 
not  be  averaged  out.  We,  however,  may  neglect  this  possibility,  because  in  the  current  work  we  consider  only 
low-inclined  moons. 

Another  situation,  which  we  exclude  from  our  treatment,  is  libration  of  the  satellite’s  periapse  about  90° 
or  270°,  caused  by  the  pull  of  a  third  body  (the  star  or  some  large  neighboring  planet).  The  possibility  of 
such  librations  may  be  derived  from  the  presence  of  the  cos  2 co  term  on  the  right-hand  side  of  the  equation 
for  da)/dt  in  the  theory  of  Kozai  (1959b,  1962) — for  an  easy  introduction  into  this  theory  see  Innanen  et  al. 
(1997),  and  for  its  generalisation  to  finite  obliquities  see  Gurfil  et  al.  (2007).  An  important  special  case  of 
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Thus,  we  arrive  at: 


^  2  R  J-,  ( l-m)\  ^2  ^2 

T  =  2^2 k,GMsal^j^  2l  ,,,  X  Flmp(0  X  Glpq(e)smelmpq  +  D  (33) 

1=2  a  m=0U  +  p= 0  q=-oo 

the  sum  standing  for  the  constant  (M -independent)  part  of  the  torque,  and  f  denoting  the 
oscillating  part  whose  time-average  is  zero. 

As  we  pointed  in  the  end  of  Sect.  5,  the  sign  of  the  phase  lag  6impq  depends  on  whether 
mO  falls  short  of  or  exceeds  the  linear  combination  (Z  —  2p)a>*  +  (Z  —  2p  +  q)M*  +  mQ.*  « 
(l  —  2p  +  q)M* .  Now  we  also  understand  that,  outside  resonances,  the  Impq  component 
of  the  tidal  torque  experienced  by  the  planet  is  decelerating  if  the  values  of  md  exceed  the 
given  combination,  and  is  accelerating  otherwise. 

Expression  (33)  gets  considerably  simplified  if  we  restrict  ourselves  to  the  case  of  l  =  2. 
Since  0  <  m  <  Z,  and  since  m  enters  the  expansion  as  a  multiplier,  we  see  that  only  m  =  1,2 
actually  matter.  As  0  <  p  <  Z,  we  are  left  with  only  six  relevant  F’s,  those  corresponding 
to  (Imp)  =  (210),  (211),  (212),  (220),  (221),  and  (222).  By  a  direct  inspection  of  the  table  of 
7qmp  we  find  that  five  of  these  six  functions  happen  to  be  0(i)  or  0(i* 2),  the  sixth  one  being 
F220  =  \  (1  +  cos  Z)2  =  3  +  0(i2 *).  Thus,  in  the  leading  order  of  i,  the  constant  part  of  the 
torque  reads: 

3  00 

Ti=2  =  2  X  GMltRa-6 *Gj0q(e)k2sme220q  +  O(i2/Q).  (34) 

q=—o o 

This  is  what  is  called  Darwin-Kaula-Goldreich  torque,  or  simply  Darwin  torque.  The  prin¬ 
cipal  term  of  this  series  is 


T2200  = 


- GM2,k2 R  a  6  sin  62200 ■ 


Switching  from  the  lags  to  quality  factors  via  formula8 


(35) 


(2lmpq  —  |C0t€impq|, 


(36) 


Footnote  7  continued 

the  theory  of  the  third-body-caused  librations  is  the  one  of  the  satellite  getting  into  a  resonance  with  the  third 

body.  An  indication  that  such  resonances  may  cause  the  satellite’s  co  librate  comes  from  the  mathematically 

similar  theory  of  Pluto-Neptune  resonances  (Williams  and  Benson  1971):  being  in  resonance  with  Neptune, 

Pluto  has  its  periapse  librating  due  to  a  high  inclination.  (To  be  exact,  the  behaviour  of  Pluto’s  periapse  is 

dictated  not  only  by  Neptune,  but  by  the  combined  influence  of  all  of  the  four  gas  giant  planets.  However,  this 
does  not  change  the  main  point:  the  outer  body  or  bodies  can  cause  apsidal  libration.) 

In  the  Solar  system,  none  of  the  large  satellites  is  so  highly  perturbed  as  to  have  a  periapse  librating  around 
90°  or  270°  due  to  the  above  two  mechanisms.  In  theory,  though,  this  remains  an  option  for  exoplanets.  Either 
librating  mechanism  might  apply  also  to  satellites  of  minor  planets.  In  our  current  paper  we  do  not  consider 
such  moons. 

8  The  phase  lag  6impq  is  introduced  in  (23-24),  while  the  tidal  harmonic  cuimpq  is  given  by  (25).  The  quality 
factor  (2lmpq  =  I  c°t^impql  is>  f°r  physical  reasons,  positively  defined.  Hence  the  multiplier  sgn  &>impq  in 
(37).  (As  ever,  the  function  sgn  (x)  is  defined  to  assume  the  values  +1,  —1,  or  0  for  positive,  negative,  or 
vanishing  x,  correspondingly.) 

Mind  that  no  factor  of  two  appears  in  (36-37),  because  €  is  a  phase  lag,  not  a  geometric  angle. 
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we  obtain: 


sin  Grnpq 


whence 


sin  |  Gmpq  I  Sgn  ^lmpq  — 


Sgn  <^>lmpq 

yj  1  T  COt^  Gmpq 


Sgn  (Dimpq  o 

&  ™pq  +  0(Q- 3), 


Q\ 


Impq 


Sgn  &>lmpq 
J1  +  Gimpq 


(37) 


3  7  5  6  2  S^n  (l)220q  2  3 

r,=2  =  2  Z  GM^R  a  G2o q(e)k2~d - “  +  0(1  /Q)  +  °(G  >• 

q=-00  ^220 q 

Now,  let  us  simplify  the  sign  multiplier.  If  in  expression  (25)  for  &>lmpq  we  get  rid  of  the 
redundant  asterisks,  replace9  M  with  Mq  +  n  ~  n,  and  set  l  =  m  =  2  and  p  =  0,  the 
outcome  will  be: 


sgn  u>210q  =  sSn  [2«  +  (2  +  <?)«  +  2Q  —  20] 

=  sgn  |^<i>  +  ^1  +  |  ]  n  +  S2  -  . 

As  the  node  and  periapse  precessions  are  slow,  the  above  expression  may  be  simplified  to 

sgn  [(l  +  f )n-e\. 

All  in  all,  the  approximation  for  the  constant  part  of  the  torque  assumes  the  form: 

3  oo  ^ 

r/=2=-  Z  GMitRa^G220 q(e)k2Q^0q 

q=—oo 

X  sgn  [(l  +  |)  n  -  0]  +  Oir/Q)  +  0{Q^).  (38) 

That  the  sign  of  the  right-hand  side  in  the  above  formula  is  correct  can  be  checked  through  the 
following  obvious  observation:  for  a  sufficiently  high  spin  rate  6  of  the  planet,  the  multiplier 
sgn  [(l  +  |)  n  ~  t>ecomes  negative.  Thereby  the  overall  expression  for  r/=2  acquires  a 
“minus”  sign,  so  that  the  torque  points  out  in  the  direction  of  rotation  opposite  to  the  direction 
of  increase  of  the  sidereal  angle  0.  This  is  exactly  how  it  should  be,  because  for  a  fixed  q 
and  a  sufficiently  fast  spin  the  q’s  component  of  the  tidal  torque  must  be  decelerating  and 
driving  the  planet  to  synchronous  rotation. 


^  While  in  the  undisturbed  two-body  setting  At  =  A4q  +  n(t  —  to)  and  At  =  n,  under  perturbation  these 
relations  get  altered.  One  possibility  is  to  introduce  (following  Tisserand  1 896)  an  osculating  mean  motion 
n(t )  =  J /i/a(t)3,  and  to  stick  to  this  definition  under  perturbation.  Then  the  mean  anomaly  will  evolve  as 
Ai  =A4o (t)  +  ft  n(t)dt ,  whence  At  =  Ato (t)  +  n(t). 

Other  possibilities  include  introducing  an  apparent  mean  motion,  i.e.,  defining  n  either  as  the  mean-anom¬ 
aly  rate  d At /dt,  or  as  the  mean-longitude  rate  dL/dt  =  dQ/dt  +  dcn/dt- 1-  d A /dt  (as  was  done  by  Williams 
et  al.  2001).  It  should  be  mentioned  in  this  regard  that,  while  the  first-order  perturbations  in  a{t)  and  in  the 
osculating  mean  motion  y/ pt/a{t )3  do  not  have  constant  parts  leading  to  secular  rates,  the  epoch  terms  typically 
do  have  secular  rates.  These  considerations  explain  why  there  exists  a  difference  between  the  apparent  mean 
motion  defined  as  dL/dt  (or  as  dAi/dt)  and  the  osculating  mean  motion  y/ [x/a{t')^ . 

In  many  practical  situations,  the  secular  rate  in  A4q  is  of  the  order  of  the  periapse  rate,  while  the  secular 
rate  in  Lo  turns  out  to  be  smaller.  Hence  the  advantage  of  defining  the  apparent  n  as  the  mean-longitude  rate 
dL/dt,  rather  than  as  the  mean-anomaly  rate  dAd/dt.  (At  the  same  time,  for  a  satellite  orbiting  an  oblate 
planet  the  secular  rates  of  Aio,  Lq,  and  periapse  are  of  the  same  order.) 

Although  the  causes  of  orbit  perturbations  are  beyond  the  scope  of  our  paper,  we  would  mention  that  in  the 
expression  (25)  for  <^jmpq  the  notations  At,  co,  and  £2  generally  imply  the  secular  rate. 
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Expansion  (38)  was  written  down  for  the  first  time,  without  proof,  by  Goldreich  and  Peale 
(1996).  A  schematic  proof  was  later  offered  by  Dobrovolskis  (2007). 


7  The  MacDonald  expression  for  the  tidal  torque 

The  idea  of  representing  the  tidal  pattern  with  one  bulge  is  often  mis-attributed  to  MacDonald 
(1964)  who  in  fact  borrowed  it  from  Gerstenkorn  (1955).  This  approach  was  furthered  by 
Singer  (1968)  and  greatly  advanced  by  Mignard  (1979, 1980).  It  is  the  latter  two  authors  who 
realised  that  Gerstenkorn’s  single-bulge  simplification  was  acceptable  only  with  a  frequency- 
independent  At,  not  with  a  frequency-independent  Q  as  in  MacDonald  (1964).  Nevertheless 
we  shall  call  this  approach  “the  MacDonald  torque”,  to  comply  with  the  established  conven¬ 
tion.  For  the  same  reason,  the  afore-described  Darwin-Kaula-Goldreich  expansion  will  be 
referred  to  simply  as  “the  Darwin  torque”. 

In  the  preceding  section,  the  Darwin  torque’s  component  orthogonal  to  the  equator  was 
conveniently  given  by  the  fundamental  formula  (29).  Within  the  MacDonald  approach,  it 
will  be  more  practical  to  write  the  torque  as  a  derivative  taken  with  respect  to  the  longitude. 
The  torque  acting  on  the  tidally  disturbed  satellite  of  mass  Msat  is  —  MsatdU /dX,  while  the 
torque  that  this  moon  exerts  on  the  planet  is  this  expression’s  negative: 

d  U(r) 

r  (r)  =  Msat— fA.  (39) 

oX 

Speaking  rigorously,  the  formula  furnishes  the  torque’s  component  perpendicular  to  the  plan¬ 
etary  equator.  As  can  be  seen  from  (44),  formula  (39)  coincides  with  (29)  for  low  inclinations. 


7.1  Simplifications  available  for  low  i 
In  principle,  we  can  as  well  insert  into 
U(r)  = 


“  /8\m 
2 ‘>(7) 

1=2  x  ' 


W,(R,r  ) 


the  “raw”  expression  (14),  the  one  as  yet  “unprocessed”  by  (15).  This  will  give  us 


U(r)  =  -GM*at^ki 


R 


1=2 


■r—i  ( l  —  m)\ 

X  7T7 - 77(2  —  <$0m) 


r”  -  (l  +  m)\ 
m= 0 


x  P/,„(sin0)P;m(sin0*)  cos/nfA.  —  X*) 
or,  for  low  inclinations  of  both  the  tidally-perturbed  and  tide-raising  satellites: 

l 


U(r)  =  —GM*al  X  k, 


R~ 


■==  ( l  —  m)\ 

X  7T7 - 77  (2  -  «Om)flm(0) 


1=2 


r,+,r*M  ,^0  (1  +  <n)\ 


x  Pim  (0)  cos  m (X  -  X*)  +  0(iz )  +  0(i*2)  +  0(ii*). 


(40) 


(41) 


At  this  point,  we  once  again  are  faced  with  the  question  of  how  to  bring  damping  into  the 
picture,  i.e.,  how  to  take  care  of  the  delayed  reaction  of  the  planet’s  material  to  the  tidal  stress. 
It  is  tempting  to  substitute  mX  with  its  delayed  value.  Then  instead  of  cos  m  (X  —  X*)  we  get 

/  (delayed)  ' 


cos  I  mX  —  mX 


^  =  cos  ^ 


mX  — 


mX  —  m  X  At 


])' 


(42) 
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This  trick,  suggested  by  Kaula  (1968,  p.  201), 10  has  a  physical  justification  only  if  At  is 
the  same  for  all  frequencies,  a  model  pioneered  by  Singer  (1968)  and  furthered  by  Mignard 
(1979,  1980).  It  can  be  shown  that  this  model  is  equivalent  to  the  following  rheological 
law:11 


Q 


lmpq 


l 

XtapqA  f  ' 


(43) 


Even  then,  though,  it  remains  unclear  how  to  connect  the  longitude  lag  m  X  *  At  with  one 
or  another  (9,,  in  the  spirit  of  (37).  To  see  what  can  be  done,  write  down  the  longitude 
(reckoned  from  a  fixed  meridian  on  the  rotating  planet)  as 

X  =  — 0  T  £2  +  o)  T  v  T  19  ( /  ' )  =  — 6  T  £2  +  co  -f-  AA 

+2e  sin  M  +  O  (e2)  +  O (i2) ,  (44) 


v  being  the  true  anomaly.  Thence,  in  neglect  of  the  nodal  and  apsidal  precessions,  the  cosine 
becomes: 

cos  —  mX  j  +  m  X  A/^  =  cos  (]^mX  —  mX  J  +  m  [v*  —  0*]  Ar^,  (45) 


or,  equivalently: 


.  ^JmX  —  mX  j  +  m  X  A =  cos  {^mX  —  mX  j  +  m  [n*  —  (?*]  At 

>) 


+  2me  n  At  cos  A)  +  0(e~ ) 


Insertion  of  (45)  into  (41),  along  with  substitution  of  r*(t)  by  r*(t  —  At),  leads  us  to 


(46) 


U(?)  =  -GM*atJ^k, 


FT 


^  r(t)'+‘r\t  -  Atj+‘  ^{l  +  m)\ 


^  (Z-w)! 

X  7^ - 7t(2-8Om)Plm(0) 


1=2 


n= 0 


x  P/,„( 0)  cos  (m  |^X  —  X  j  +  m  [v*  —  0*]  Ar^  +  0(i2)  +  0(i*2)  +  0(ii*). 


(47) 


If  we  take  into  account  only  the  l  =  2  contribution,  expression  (41)  will  simplify  to 


GM*k2R ‘  ^  (2  —  m)\ 

C/(r)  = - X  7^ - 77  (2  -  8om)P2m(0)P2m(0) 


r(t)  r  (?) 


12=0 


(2  +  m)\ 


x  cosra(A.  —  A*)  +  0(il)  +  0(i*2)  +  0(ii*) 


(48) 


10  Mind  the  difference  in  notations.  While  in  the  original  paper  Kaula  (1964)  denoted  the  phase  lags  with 
eimpq » his  book  Kaula  (1968)  called  them  ^lmpq .  For  the  longitudinal  lag  2  X  A t  emerging  in  our  formula 
(46),  Kaula  (1968)  used  notation  28.  This  way,  in  the  terms  used  by  Kaula  (1968)  in  his  book,  the  geometric 
angle  subtended  at  the  primary’s  center  between  the  directions  to  the  bulge  and  the  moon  is  called  28,  not  8 
as  in  most  literature. 

11  Combining  (24)  with  the  relation  Q  =  1/tan  |e|,  we  see  that  setting  all  A/lmpq  equal  to  the  same  At  is 
equivalent  to  saying  that  the  quality  factor  scales  as  the  inverse  frequency:  Q  =  1  / (x  AO,  provided,  of  course, 
that  the  Q  factor  is  large.  As  can  be  seen  from  (37),  a  more  exact  relation  will  read:  sin(x  AO  =  l/\/l  +  <22, 
so  that  x  A t  =  Q~l  +  0{Q~ 3). 

Very  special  is  the  case  when  the  values  of  the  quality  factor  are  very  low  (say,  much  less  than  10).  In  this 
situation,  the  interconnection  between  the  quality  factor  and  the  phase  lag  becomes  quite  different  from  the 
customary  formula  Q  =  cot  |e|.  See  the  Appendix  for  details. 
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where  only  the  m  =  2  term  is  important.12  In  the  presence  of  dissipation,  the  appropriately 
simplified  version  of  (48)  will  read: 


U(r)  =  -- 


3  GM*k2R 


4  r(t)  r*(t  -  At) 


:  (J2A  -  2A.*]  +  2  [v*  -  9*]  A t) 


+  0(i2)  +  0(i*2)  +  0(ii*), 


(49) 


while  the  corresponding  expression  for  the  torque  exerted  by  the  satellite  on  the  planet  will, 
in  this  approximation,  be  given  by 


r  (r)  =  MSl 


9 U(r)  3  GM*Msatk2R 


dX 


2  r(t)  r*(t  —  Aty 


sin  ^2X  -  2A*]  +  2  [v*  -  9*]  A t) 


+  0(i2)  +  O  (/*“)  +  0(ii*). 


(50) 


In  the  case  when  the  tidally  disturbed  satellite  coincides  with  the  tide-raising  one,  i.e.,  when 
X  =  X*  and  Msat  =  M*at,  we  obtain: 

3  R 5 

r  =  -GM2.dtk2—— - — T  sin  (2  [v  -  9]  At)  +  0(i2/Q) 

Z  r{t)  r(t  —  At) 

3  R 5 

=  -GM^k 2—  sin  (2  [i.  -  9]  At)  +  0(?/Q)  +  0(en/Q2X),  (51) 

z  r 

where  the  error  Oien/ Q2X)  emerges  when  we  identify  the  lagging  distance  r(t  —  At)  with 
r  =  r(t).  Replacement  of  r{t  —  At)  with  r  is  convenient,  though  not  necessary.  In  Subsect. 
7.2  below,  we  shall  explain  that,  after  averaging  over  one  revolution  of  the  moon  about  the 
planet,  the  error  caused  by  this  replacement  reduces  to  O (e2 n2 / Q3 X2) ,  which  will  be  less 
than  the  largest  error. 

The  MacDonald  torque  (51)  is  equivalent  to  the  Darwin  torque  (34)  with  an  important 
proviso  that  all  time  lags  Aflmpq  are  equal  to  one  another  or,  equivalently,  that  the  rheological 
model  (43)  is  accepted.  Physically,  the  special  case  of  equal  time  lags  is  exactly  the  case  when 
the  tide  may  be  rigorously  interpreted  as  one  double  bulge  of  a  variable  rate  and  amplitude. l  j 
Mathematically,  this  model  enables  one  to  wrap  up  the  infinite  series  (34)  into  the  elegant 
finite  form  (51).  Formally,  this  wrapping  can  be  described  like  this:  expression  (51)  mimics 
the  principal  term  of  the  series  (34),  provided  in  this  term  the  multiplier  G700  is  replaced 
with  unity,  a  is  replaced  with  r,  and  the  principal  phase  lag 

622oo  “  2(”  —  9)  At  (52) 

is  replaced  with  the  longitudinal  lag  or,  possibly  better  to  say,  with  the  quasi-phase 

€  =  2(i)  —  9)  At.  (53) 


Thus,  we  see  that  within  the  MacDonald  one-variable-bulge  formalism  the  longitudinal  lag 
(53)  is  acting  as  an  instantaneous  phase  lag  associated  with  the  instantaneous  tidal  frequency 
X  =  2|v  —  9 1.  This  is  why  we  may  call  it  simply  6,  without  a  subscript.  Evidently,  e  is  (up  to 


12  In  (48),  we  may  neglect  the  ^-independent  term  with  m  =  0,  because  our  eventual  intention  is  to  find  the 
torque  by  differentiating  U (r)  with  respect  to  X.  We  may  also  omit  the  m  =  1  term,  because  P2\ (0)  =  0.  This 
omission  brings  up  an  error  of  order  0(ii*)  into  equations  (41),  (47-50)  and  (56). 

12  An  attempt  to  generalise  this  simplified  approach  to  arbitrary  inclinations  was  undertaken  by  Efroimsky 
(2006).  While  for  constant  time  lags  that  generalisation  is  likely  to  be  acceptable,  it  remains  to  be  explored 
whether  the  offers  a  practical  approximation  for  actual  rheologies  (60). 


Springer 


272 


M.  Efroimsky,  J.  G.  Williams 


a  sign)  twice  the  geometrical  angle  subtended  at  the  primary’s  center  between  the  directions 
to  the  moon  and  to  the  bulge. 14 

The  geometric  meaning  of  the  longitudinal  lag  being  clear,  let  us  consider  its  physical 
meaning,  in  the  sense  of  this  lag’s  relation  to  the  dissipation  rate.  For  some  fixed  frequency 
Ximpq,  the  corresponding  phase  lag  6impq  is  related  to  the  appropriate  quality  factor  via 
1/gimpq  =  tan  |eimpq|.  To  keep  the  analogy  between  the  true  lags  and  the  instantaneous  lag 
(53),  one  may  conveniently  define  a  quantity  Q  as  the  inverse  of  tan  |e|.  This  will  enable  one 
to  express  the  MacDonald  torque  as 

3  R 5 

r  =  ^GMs2atk2— 3— -  3  sine  +  O^/Q) 

^  r(t)  r(t  —  At) 

3 

=  -GM2atk2  —  Q-lsgn(v-9)+0(i2/Q)  +  0(en/Q2X)+0(Q-3)  (54) 

2  r 

Since  Q  was  defined  as  l/tan|e|,  it  is  not  guaranteed  to  deserve  the  name  of  an  overall 
quality  factor.  At  each  particular  frequency  Ximpq ,  the  corresponding  quality  factor  gimpq  = 
1/tan  |eimpq]  is  related  to  the  peak  energy  of  this  mode,  TipeakCXimpq)’  and  to  the  one-cycle 
energy  loss  at  this  frequency,  A£cyc|e(/lmpq),  via 


A  £Cycle  (  Xlinpq  ) 


£peak  (Ximpq  ) 


Q 


lmpq 


(55) 


However,  it  is  not  at  all  obvious  if  the  quantity  Q  defined  through  the  longitudinal  lag  as 
Q  =  1/tan  |e|  interconnects  the  overall  tidal  energy  with  the  overall  one-cycle  loss,  in  a 
manner  similar  to  (55).  The  literature  hitherto  has  always  taken  for  granted  that  it  does. 
However,  the  proof  (to  be  presented  elsewhere)  requires  some  effort.  The  proof  is  based  on 
interpreting  y  =  2|v  —  6\  as  an  instantaneous  tidal  frequency. 

The  interconnection  between  Q  =  1/  tan  \f  j  and  the  overall  energy-damping  rate  mimics 
(55)  only  up  to  a  relative  error  of  order  0(en/ Qy)  =  O(enAt),  i.e.,  up  to  an  absolute  error 
of  order  0(eQ~1nAt).  This  is  acceptable,  because  in  realistic  settings  n  At  1. 


7.2  Further  simplifications  available  in  the  zeroth  order  of  en/ Qy 

Suppose,  we  ignore  the  difference  between  r  and  r*,  which  are  the  two  locations  of  the  same 
satellite,  separated  by  the  time  lag  owing  to  the  tidal  response.  We  shall  now  demonstrate  that, 
though  the  relative  error  of  this  approximations  is  0(en/  Qy),  after  averaging  over  a  satellite 
period  this  approximation  brings  only  a  0(e2n2/Q2y2)  relative  error  into  the  expression  for 
the  torque. 

From  the  well-known  formulae  r  =  a(  1  —  e2)/(l  +  ecosv)  and  3 v/dM  =  (1  + 
ecosv)2/(l  —  e2)3/2  we  see  that 


A r  =  r(t )  -  r(t  -  At)  =  - 


ae(  1  —  e2) 
(1  +  ecos  v)2 


sin  v  An  -|-  O  (e(  Av)2) 


ae  sin  v 

( 1  —  e2)1/2 


n  At  +  O  (e(nAt)2)  . 


14  As  the  subtended  angle  is  |(v  —  3)A/|.  its  double  is  equal  to  the  absolute  value  of  e,  and  not  to  that  of 
<2204  =2("-^ 
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The  time  lag  is  interconnected  with  the  phase  shift  and  the  quality  factor  via  the  relations 

XAt  =  e*  Q~\ 

X  =  2\0  —  v\  being  the  instantaneous  tidal  frequency.  Hence 

c  n 

A r  =  r(t)  —  r(t  —  At)  «  —a - sin  v. 

Q  x 

As  A r  is  proportional  to  sin  v,  only  terms  quadratic  in  A r  survive  averaging.  Thus,  while  in 
U  =  cos(2A.  -  2>*  +  e)  +  0(en/QX)  +  0(i2)  +  0(i*2)  +  0(ii*), 


(56) 


and 


(57) 


the  relative  error  is  0(en/QX )  +  0(i2),  in  the  averaged  expression15 


sin  e)+ 0(e2n2/Q3X2)+ 0(i2/Q) 


(58a) 


3  GM2satk2R 


■2n  R4 

—t  sine  dv  +  O (e2 n2 / Q3 X2)  +  0(i2/Q).  (58b) 
r4 


1 


L 


Ana2 


(l-e2)1/2 


it  is  only  0(e2n2 / Q2X2)  +  0(i2). 

In  the  above  expressions,  we  asserted  after  the  differentiation  that  M*at  =  Msa t  and 
K*  =  X,  implying  that  the  tide-generating  and  tidally-perturbed  moons  are  one  and  the  same 
body.  As  soon  as  X  is  set  to  be  equal  to  X*,  the  sine  function  in  (58)  becomes  sine  ~  l/<2- 
So,  while  the  relative  error  in  (58)  is  0{e2n2 / Q2x2)  +  0(i2),  the  absolute  error  becomes 
0(eV/Q3X2)  +  Oy/Q). 

The  error  0(e2n2 / Q3 X2)  becomes  irrelevant  for  two  reasons.  First,  our  substitution  of 
sine  with  tane  =  1  / Q  generates  a  relative  error  0(Q  2 ) ,  i.e.,  an  absolute  error  0(Q  3 ) . 


15  We  recall  that  time  averages  over  one  revolution  of  the  satellite  about  the  primary  are  given  by 


while  the  planetocentric  distance  isr=a(l  —  e2\  /  (\  +  e  cos  v),  with  v  being  the  true  anomaly.  This  way. 
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Second,  as  explained  in  the  end  of  Subsect.  7.1,  the  uncertainties  inherent  in  our  definition 
of  the  overall  quality  factor  Q  entail  an  absolute  error  0(en/Q2/  ).  Each  of  these  two  errors 
exceeds  0(e2n2 / Q2 x2).  We  can  then  write: 

3 GM~,ki  I  sgn  (6  —  v)  R6  \  ,  7  , 

<r>  = - \  — — ~Q - ^  +°(e'3)+O(i2/0+OWe2/)  (59a) 

3 GM2satk2R  1  r2*  R4  sgn  (0  -  v) 

Ana2  (i_  *2)1/2  ./o  74  q  dv 

+  0{Q-2)  +  0(i2/Q)  +  0(en/Q2X).  (59b) 


8  Use  and  abuse  of  approximation  (56-59) 


Just  as  with  the  formula  (56)  for  the  potential,  the  elegant  expression  (57)  for  the  torque 
remains  correct  only  to  the  zeroth  order  in  e/<2,  while  (59)  is  valid  to  the  first  order.  This 
is  the  reason  why  the  convenience  of  this  approximation  and  of  its  corollaria  is  somewhat 
deceptive.  Nevertheless,  the  (56)-(59)  were  employed  by  many  an  author. 

Goldreich  and  Peale  (1996)  used  them  to  build  a  theory  containing  terms  up  to  e1 .  Now 
we  see  that  some  coefficients  in  their  theory  of  capture  into  resonances  must  be  reconsidered. 
The  same  pertains  to  some  coefficients  in  the  theory  of  Mercury’s  rotation,  recently  offered 
by  Peale  (2005).  Fortunately,  the  key  conclusions  of  Peale  (2005)  stay  unaltered,  despite  the 
corrections  needed  in  the  said  coefficients.16 

Interestingly,  Kaula  (1968)  fell  into  this  temptation,  and  so  did  Goldreich  (1966b).  Equa¬ 
tion^. 5.29)  in  Kaula  (1968),  as  well  as  Eq.  (15)  in  Goldreich  (1966b),  is  but  the  above 
formula  (59)  with  the  inverse  quality  factor  taken  out  of  the  integral: 


^ Kaula 


3  GM2satk2R  1 

An  a2Q  (\-e2)X/2 


v)dv. 


(Kaula  1968,  Eq.  4.5.29) 


Besides  the  afore-mentioned  fact  that  this  approach  contains  a  relative  error  0(en/  Q/  )  + 
0(QT2)  +  0(i2),  it  suffers  a  greater  defect.  Taking  Q  out  of  the  integral  is  illegitimate, 
because  it  implies  frequency-independence  of  Q.  This  is  then  incompatible  with  Kaula’s 
implicit  assumption  of  a  constant  At,  an  assumption  tacitly  present  in  (42). 17 

Goldreich  (1966b)  and  Kaula  (1968)  used  this  oversimplified  formula  to  investigate  libra- 
tions  of  a  satellite  trapped  in  a  1:1  resonance.  Other  authors  used  it  to  evaluate  despinning 
rates  of  bodies  outside  this  resonance.  We  shall  dwell  on  the  latter  case  in  Sect.  10. 


9  Can  the  quality  factor  scale  as  a  positive  power  of  the  tidal  frequency? 

As  of  now,  the  functional  form  of  the  dependence  Q  (/ )  for  Jovian  planets  remains  unknown. 
For  terrestrial  planets,  the  model  Q  ~  I  //  is  definitely  incompatible  with  the  geophysical 
data.  A  convincing  volume  of  measurements  firmly  witnesses  that  Q  of  the  mantle  scales  as 
the  tidal  frequency  to  a  positive  fractional  power: 

16  Stan  Peale,  private  communication,  2007. 

17  This  would  be  incompatible  with  the  MacDonald  (1964)  treatment  as  well,  because  MacDonald’s  formal¬ 
ism  necessitates  the  rheology  Q  ~  l/x ■  We  shall  return  to  this  point  in  Subsect.  10.1. 
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Q  =  (£x)“  .  where  a  =  0.3  ±  0.1, 


(60) 


£  being  an  integral  rheological  parameter  with  dimensions  of  time.  This  rheology  is  incom¬ 
patible  with  the  postulate  of  frequency-independent  time-delay.  Therefore,  insertion  of  the 
realistic  model  (60)  into  the  formula  presented  in  Sect.  8  will  remain  insufficient.  An  honest 
calculation  should  be  based  on  averaging  the  Darwin-Kaula-Goldreich  formula  (38),  with 
the  actual  scaling  law  (60)  inserted  therein,  and  with  the  appropriate  dependence  Aflmpq  (xlmpq) 
taken  into  account  (see  formula  (83)  below). 

9. 1  The  “paradox” 

Although  among  geophysicists  the  scaling  law  (60)  has  long  become  common  knowledge, 
in  the  astronomical  community  it  is  often  met  with  prejudice.  The  prejudice  stems  from  the 
fact  that,  in  the  expression  for  the  torque,  Q  stands  in  the  denominator: 


At  the  instant  of  crossing  the  synchronous  orbit,  the  principal  tidal  frequency  x2200  becomes 
nil,  for  which  reason  insertion  of 


Q~Xa,  «>0 

into  (61)  seems  to  entail  an  infinitely  large  torque  at  the  instant  of  crossing: 


(62) 


(63) 


a  clearly  unphysical  result. 

Another,  very  similar  objection  to  (60)  originates  from  the  fact  that  the  quality  factor  is 
inversely  proportional  to  the  phase  shift:  Q  ~  1/e.  As  the  shift  (24)  vanishes  on  crossing 
the  synchronous  orbit,  one  may  think  that  the  value  of  the  quality  factor  must,  effectively, 
approach  infinity.  On  the  other  hand,  the  principal  tidal  frequency  vanishes  on  crossing  the 
synchronous  orbit,  for  which  reason  (60)  makes  the  quality  factor  vanish.  Thus,  we  come  to 
a  contradiction. 

For  these  reasons,  the  long-entrenched  opinion  is  that  these  models  introduce  discontinu¬ 
ities  into  the  equations  and  can  thus  be  considered  as  unrealistic  approximations  for  rotating 
bodies. 

It  is  indeed  true  that,  while  law  (60)  works  over  scales  shorter  than  the  Maxwell  time 
(about  102  years  for  most  minerals),  it  remains  subject  to  discussion  in  regard  to  longer 
timescales.  Nonetheless,  it  should  be  clearly  emphasised  that  the  infinities  emerging  at  the 
synchronous-orbit  crossing  can  in  no  way  disprove  any  kind  of  rheological  model.  They  can 
only  disprove  the  flawed  mathematics  whence  they  provence. 

9.2  A  case  for  reasonable  doubt 

To  evaluate  the  physical  merit  of  the  alleged  infinite-torque  “paradox”,  recall  the  definition  of 
the  quality  factor.  As  part  and  parcel  of  the  linearity  approximation,  the  overall  damping  inside 
a  body  is  expanded  in  a  sum  of  attenuation  rates  corresponding  to  each  periodic  disturbance: 


(64) 


where,  at  each  frequency 
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<£(x,)> 


(£(X,))  =  _  £peak(X,) 

Q(Xt)  Xi  G(x,) 


(65) 


(. . .)  designating  an  average  over  a  flexure  cycle,  E(y.)  denoting  the  energy  of  defor¬ 
mation  at  the  frequency  x, ,  and  <2(X,)  being  the  quality  factor  of  the  medium  at  this 
frequency. 

This  definition  by  itself  leaves  enough  room  for  doubt  in  the  above  “paradox”.  As  can 
be  seen  from  (65),  the  dissipation  rate  is  proportional  not  to  1/Q(x)  but  to  x/G(x)-  This 
way,  for  the  dependence  Q  ~  the  dissipation  rate  ( E )  will  behave  as  x1-“-  In  the  limit 
of  x  — >  0,  this  scaling  law  portends  no  visible  difficulties,  at  least  for  the  values  of  a  up  to 
unity.  While  raising  a  above  unity  may  indeed  be  problematic,  there  seem  to  be  no  funda¬ 
mental  obstacle  to  having  materials  with  positive  a  taking  values  up  to  unity.  So  far,  such 
values  of  a  have  caused  no  paradoxes,  and  there  seems  to  be  no  reason  for  any  infinities  to 
show  up. 


9.3  The  phase  shift  and  the  quality  factor 

As  another  preparatory  step,  we  recall  that,  rigorously  speaking,  the  torque  is  proportional 
not  to  the  phase  shift  e  itself  but  to  sin  e.  From  (37)  and  (60),  we  obtain: 

|  sin  el  =  —  =  =  (66) 

vi  +  q-  yi+^x- 

We  see  that  only  for  large  values  of  Q  one  can  approximate  |  sine|  with  l/Q  (crossing 
of  the  synchronous  orbit  not  being  the  case).  Generally,  in  any  expression  for  the  torque, 
the  factor  l/Q  must  always  be  replaced  with  1  fy/ 1  +  Q2.  Thus,  instead  of  (61)  we  must 
write: 

1  1 

r  ~  |  sin  6 1  =  =  =  ,  =,  (67) 

vi  +  q  yi+f2v“ 

£  being  a  dimensional  constant  from  (60). 

Although  this  immediately  spares  us  from  the  fake  infinities  at  X  — >■  0,  we  still  are  facing 
a  strange  situation:  it  follows  from  (66)  that,  for  a  positive  a  and  vanishing  the  phase 
lag  6  must  be  approaching  it /I,  thereby  inflating  the  torque  to  its  maximal  value  (while  on 
physical  grounds  the  torque  should  vanish  for  zero  x)-  Evidently,  some  important  details  are 
still  missing  from  the  picture. 

9.4  The  stone  rejected  by  the  builders 

To  find  the  missing  link,  recall  that  Kaula  (1964)  described  tidal  damping  by  employing  the 
method  suggested  by  Darwin  (1880):  he  accounted  for  attenuation  by  merely  adding  a  phase 
shift  to  every  harmonic  involved — an  empirical  approach  intended  to  make  up  for  the  lack  of 
a  consistent  hydrodynamical  treatment  with  viscosity  included.  It  should  be  said,  however, 
that,  prior  to  the  work  of  1880,  Darwin  had  published  a  less  known  article  (Darwin  1879), 
in  which  he  attempted  to  construct  a  self-consistent  theory,  one  based  on  the  viscosity  factor 
of  the  mantle,  and  not  on  empirical  phase  shifts  inserted  by  hand.  Darwin’s  conclusions  of 
1879  were  summarised  and  explained  in  a  more  general  mathematical  setting  by  Alexander 
(1973). 
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The  pivotal  result  of  the  self-consistent  hydrodynamical  study  is  the  following.  When  a 
variation  of  the  potential  of  a  tidally  distorted  planet,  U(r),  is  expanded  over  the  Legendre 
functions  P/„,(sin0),  each  term  of  this  expansion  will  acquire  not  only  a  phase  lag  but  also 
a  factor  describing  a  change  in  amplitude.  This  forgotten  factor,  derived  by  Darwin  (1879), 
is  nothing  else  but  cose.  Its  emergence  should  in  no  way  be  surprising  if  we  recall  that  the 
damped,  forced  harmonic  oscillator 

x  +  2  yx  +  o?0x  =  Fe,Xt  (68) 


evolves  as 

x(t)  =  C\e 
where  the  phase  lag  is 


(-y-H'V  r— 

+  C2e 


F  cos  e  i(rt-f) 

“b  2  ’ 

-  k~ 


tan  e  = 


2yk 


K  - x2) 


(69) 


(70) 


and  the  first  two  terms  in  (69)  are  damped  away  in  time.18 

In  the  works  by  Darwin’s  successors,  the  allegedly  irrelevant  factor  of  cos  e  fell  through 
the  cracks,  because  the  lag  was  always  asserted  to  be  small.  In  reality,  though,  each  term  in 
the  Fourier  expansions  (21),  (30-35),  and  (38)  should  be  amended  with  cos  e  .  Likewise, 

*  lmpq 

the  correct  versions  of  (50-51)  and  (54)  should  contain  an  extra  factor  of  cose  .  For  the 
v  '  2200 
same  reason,  instead  of  (67),  we  should  write  down: 


r  ~  |  cos  e  sin  e  | 


Q  _ l_ 

Vl  +  Q2  Vl  +  Q2 


•,  .  f-,2a  2a 

1  +£  x 


(71) 


At  this  point,  it  would  be  tempting  to  conclude  that,  since  (71)  vanishes  in  the  limit  of  y  — >•  0, 
for  any  sign  of  a,  then  no  paradoxes  happen  on  the  satellite’s  crossing  the  synchronous  orbit. 
Sadly,  this  straightforward  logic  would  be  too  simplistic. 

In  fact,  prior  to  saying  that  cos  e  sin  e  —>■  0,  we  must  take  into  consideration  one  more 
subtlety  missed  so  far.  As  demonstrated  in  the  Appendix,  taking  the  limit  of  Q  — »■  0  is  a 
nontrivial  procedure,  because  at  small  values  of  Q  the  interconnection  between  the  lag  and  the 
Q  factor  becomes  very  different  from  the  conventional  Q  =  cot  |e|.  A  laborious  calculation 
shows  that,  for  Q  <  1  —  ;r/4,  the  relation  becomes: 


Isinecosel  =  (3g)^3 


1-^(3  g)2/3  + 0(G4/3) 


which  indeed  vanishes  for  g  — >>  0.  Both  62200  and  the  appropriate  component  of  the  torque 
change  their  sign  on  the  satellite  crossing  the  synchronous  orbit. 

So  the  main  conclusion  remains  in  force:  nothing  wrong  happens  on  crossing  the 
synchronous  orbit,  Q.E.D. 


18  As  demonstrated  by  Alexander  (1973),  this  example  indeed  has  relevance  to  the  hydrodynamical  theory  of 
Darwin,  and  is  not  a  mere  illustration.  Alexander  (1973)  also  explained  that  the  emergence  of  the  cos  e  factor 
is  generic.  (Darwin  (1879)  had  obtained  it  in  the  simple  case  of  /  =  2  and  for  a  special  value  of  the  Love 
number:  &2  =  1-5.) 

A  further  investigation  of  this  issue  was  undertaken  in  a  comprehensive  work  by  Churkin  (1998),  which 
unfortunately  has  never  been  published  in  English  because  of  a  tragic  death  of  its  Author.  In  this  preprint, 
Churkin  explored  the  frequency-dependence  of  both  the  Love  number  &2  and  the  quality  factor  within  a  broad 
variety  of  rheological  models,  including  those  of  Maxwell  and  Voight.  It  follows  from  Churkin’s  formulae  that 
within  the  Voight  model  the  dynamical  k 2  relates  to  the  static  one  as  cos  €.  In  the  Maxwell  and  other  models, 
the  ratio  approaches  cos  e  in  the  low-frequency  limit. 
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10  Tidal  despinning 

The  following  formula  for  the  average  deceleration  rate  6  of  a  planet  due  to  a  tide-raising 
satellite  has  often  appeared  in  the  literature: 

(0)  =  -1C  [9A{e)  -  nU{e)]  ,  (72) 

where 


and 

Nie)  =  (l  +  ^e2  +  ye4  +  yU6)  (l  -  e2f6  ,  (74) 

6  being  the  sidereal  angle,  8  being  the  primary’s  spin  rate,  K,  being  some  constant,  and 
the  angular  brackets  designating  an  average  over  one  revolution  of  the  secondary  about  the 
primary.  This  expression  was  derived  by  different  methods  in  Goldreich  and  Peale  (1996)  and 
Hut  (1981),  and  was  later  employed  by  Dobrovolskis  (1995,  2007)  and  Correia  and  Laskar 
(2004). 19 

10.1  Derivation  by  means  of  the  MacDonald  torque 

The  following  proof  of  (72-74)  is  implied  in  Goldreich  and  Peale  (1996)  and  is  presented  in 
more  detail  in  Dobrovolskis  (2007).  Their  starting  point  was  the  MacDonald  torque  (59).  Hut 
(1981),  who  approached  the  issue  in  the  language  of  the  Lagrange-type  planetary  equations, 
took  into  account,  in  the  disturbing  function,  only  the  leading  term  of  series  (21),  and  substi¬ 
tuted  the  principal  tidal  frequency  x2200  =  2|0  —  n |  with  the  synodal  frequency  x  =  2\G  —  v\. 
Thereby,  his  approach  was  equivalent  to  that  of  Goldreich  and  Peale  (1996)  and  Dobrovolskis 
(2007). 

Although  not  necessarily  assumed  by  these  authors,20  their  method,  as  we  saw  in  the 
section  above,  inherently  implied  the  following  assertions: 

(I)  The  quantity  x  =  2|0  —  v|  is  treated  as  an  instantaneous  tidal  frequency.  Accordingly, 
the  overall  quality  factor  Q  is  implied  to  be  a  function  not  of  the  principal  frequency 
X2200  but  of  the  instantaneous  frequency  x  ■ 

(II)  The  functional  form  of  this  dependence  is  chosen  as  Q  =  (At)-1/^1,  where  At  is 
the  time  lag. 

(Ill)  The  time  lag  At  is  frequency-independent.  This  assertion  is  equivalent  to  (II),  as  can 
be  demonstrated  from  (24). 


19  Our  formula  (78)  differs  from  formula  (4)  in  Correia  and  Laskar  (2004)  by  a  factor  of  n/x,  because  in  Ibid. 
the  quality  factor  had  been  introduced  as  1  /  (n  A  t )  and  notas  1 / (/ A  r )  ■ 

20  It  should  be  mentioned  that  the  original  treatment  by  MacDonald  (1964)  is  inherently  contradictory.  On 
the  one  hand,  MacDonald  postulates  (following  Gerstenkorn  1955)  that  there  exists  one  overall  double  bulge. 
As  explained  in  Subsect.  7.1  above,  this  assertion  unavoidably  implies  constancy  of  the  time  lag  At,  so  that 
Q  ~  l/X  and  e  ~  x-  However,  MacDonald  (1964)  erroneously  set  Q  (and,  thence,  also  f  )  frequency-inde¬ 
pendent,  an  assertion  incompatible  with  his  and  Gerstenkorn’s  postulate  of  existence  of  an  overall  double 
bulge. 

Whenever  in  the  current  paper  we  refer  to  MacDonald’s  torque,  we  always  imply  his  postulate  that  one 
double  bulge  exists.  At  the  same  time,  to  make  the  MacDonald-Gerstenkom  treatment  consistent,  we  always 
adjust  the  MacDonald-Gerstenkom  treatment  by  letting  Q  and  e  scale  as  1  /y  and  / ,  correspondingly. 


(I-,2) 


-9/2 


(73) 


—  j  1  -f-  3e~  + 
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Beside  this,  those  authors  neglected  the  order -en/ Qx  difference  between  r  and  r*  in  (54), 
generating  a  relative  error  in  r  of  order  0(en/ Qx)  (which,  luckily,  reduced  to  0(e2n2/ Q2x2) 
after  orbital  averaging).  They  also  substituted  sin  e  with  1  /  Q ,  causing  a  relative  error  of  order 
0(1/  Q2),  because  in  reality  Q  is  the  reciprocal  of  tan  e,  not  of  sine. 

Assertion  (II)  can  be  written  down  in  more  generic  notation: 

Q  =  (£x)  ,  with  a  =  —1.  (75) 


This  form  of  the  scaling  law  is  more  convenient,  for  it  leaves  one  an  opportunity  to  switch 
to  different  values  of  a.  For  any  value  of  a  (not  only  for  —1),  the  constant  £  is  an  integral 
rheological  parameter  (with  the  dimension  of  time),  whose  physical  meaning  is  explained 
in  Efroimsky  and  Lainey  (2007).  It  can  be  shown  that  in  the  particular  case  of  a  =  —  1  the 
parameter  £  coincides  with  At.  In  realistic  situations,  a  differs  from  —  1,  while  the  parameter 
£  is  related  to  the  time  lag  in  a  more  sophisticated  way  (Ibid.). 

To  show  how  (72-74)  stem  from  the  above  Assertions,  keep  for  the  time  being  a  =  —  1. 
Also  recall  that  the  torque  is  despinning  (so  9  >  n),  and  that  for  the  averages  over  time 

(r) 

(0)  =  ^r.  (76) 

C  =  f  Mpianet^2  being  the  maximal  moment  of  inertia  of  the  planet.  (For  a  homogeneous 
spherical  planet,  £  =  2/5.)  Then  plug  (75)  into  (59)  and  average  the  torque:21 


<r>  =  - 


3  GMSSLlk2£ 

h 


R 


(9  -  vW  )  +  0(iL/Q)  +  0(Q ~3)  +  O (en / Qz x) 


3  GMltk2£jtf\  +  3  GMlMk2£  f.R 


R 


R 


+  0(i~/Q)  +  0(Q->)  +  0(en/Q~x) 

3 GM2satk2£  Rb  -  i  ^ 

^6 


R 


ou 


9/2  1 

2:r  jo 

2  71 


(1  +  ecos  vy  dv 


+  3GMltk2£tf  2  -6  J_  , 

R  a6  '  In  J o 

+  0(i2/Q)  +  0(Q-3)  +  O  (en  /  Q2  X) , 


(1  +  ecos  vy  dv 


(77a) 


(77b) 


where  the  absolute  error  O (en /  Q2 x)  emerges  due  to  an  uncertainty  in  the  definition  of  the 
overall  quality  factor  Q  employed  in  MacDonald’s  model. 

Evaluation  of  the  above  integrals  is  trivial  and  indeed  leads  to  (12-1  A),  the  constant  being 

3G M/tak.2£  R6  ?>n2 M/^kiAt  R 

CR  pianet  (Mpianet  +  ;,17sat)  a3 

=  _ 3»Msatfc2 _ R^n 

£  Q 37p]anet ( T/planer  4“  4Tsa t)  a3  X 

where  we  used  the  fact  that  for  a  =  —  1  the  rheological  parameter  £  is  simply  the  lag  At. 


21  As  explained  in  the  paragraph  preceding  formula  (59),  substitution  of  sin  e  with  I  /  Q  i  n  the  expression  for 
torque  generates  a  relative  error  0(Q~2),  i.e.,  an  absolute  error  0(Q~3).  Instead  of  inserting  (75)  into  (59), 
one  may  directly  use  (53).  Still,  approximation  of  sin  e  with  e  will  entail,  in  (77)  and  its  corollaria,  a  relative 
error  0(Q~2)  and  an  absolute  error  0(Q~3).  The  situation  will  become  more  complicated  in  the  special  case 
of  low  values  of  the  quality  factor.  See  the  Appendix  below. 
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It  should  also  be  added  that,  since  (77b)  contains  a  relative  error  0(Q  2),  the  usefulness 


of  the  e4  and  eb  terms  in  (73-74)  depends  on  the  values  of  the  eccentricity  and  the  quality 


factor.  If,  for  example,  Q  =  70,  then  the  e4  terms  become  unimportant  for  e  <  0.12,  while 
the  e 6  terms  become  unimportant  for  e  <  0.24. 

To  draw  to  a  close,  we  would  mention  that  besides  the  above  formula  (72),  in  the  litera¬ 
ture  hitherto  we  saw  its  sibling,  an  expression  derived  in  a  similar  way,  but  with  Assertion 
II  rejected  in  favor  of  treating  Q  as  a  frequency-independent  constant.  The  result  of  this 
treatment  suffers  an  incurable  birth  trauma — the  incompatibility  between  the  frequency- 
independence  of  At  and  the  frequency-independence  of  Q. 

10.2  Calculation  based  on  the  Darwin  torque 

The  following  alternative  derivation  is  based  on  the  same  Assertions  (I— III)  and,  naturally, 
leads  to  the  same  results.  The  idea  is  to  calculate  the  despinning  rate  not  in  terms  of  the  Mac¬ 
Donald  torque,  but  in  terms  of  the  Darwin  torque,  keeping  the  eccentricity-caused  relative 
error  at  the  level  of  0(eb). 

To  keep  the  inclination-caused  relative  error  at  the  level  of  0{i2),  we  still  assume,  in  (34), 
that  /  =  2,  m  =  2,  p  =  0.  As  for  the  the  values  of  q,  we  keep  only  the  ones  giving  us 
terms  of  order  up  to  e  ,  inclusively.  Besides,  we  assume  the  phase  lags  to  be  small,  so  that 
sine^  =  6lmpq  +  0(e3)  =  flmpq  +  0(<2~3).  Under  all  these  presumptions,  the  constant  part 
of  the  tidal  torque  can  be  approximated  with 


2 


3 

2 


GMltR5a~6k2  Y  G~  S3n e22oq  +  °(e6/2)  +  0(/2/G)  (79a) 


r 


1=2 


q=- 2 
2 


3 

2 


GMltR5a~6k 2  Y  G]oqe220 q  +  °(e6/2)  +  OU2/Q)  +  0(Q~ 3),  (79b) 


where,  according  to  the  tables  (Kaula  1966), 


(80) 


and,  according  to  formula  (24), 


e. 


220-2 


6, 


(81) 


e 


Provided  the  quality  factor  scales  as  inverse  frequency,  all  the  time  lags  are  the  same  con¬ 
stant  At,  so  the  above  formulae  all  together  entail,  in  the  case  of  nonresonant  prograde 
spin: 


(82) 
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which  coincides  with  (72-74)  to  the  order  e 4,  inclusively,  provided  we  substitute  X2200  instead 
of  x  in  the  expression  (78)  for  K,. 


10.3  Rheologies  different  from  Q  ~  1/x 


A  part  and  parcel  of  both  afore-presented  methods  was  the  assertion  of  all  the  time  lags  Arlmpq 
being  equal.  In  reality,  the  time  lags  vary  from  one  harmonic  to  another. 

Any  particular  functional  form  of  the  dependence  At  (/)  fixes  the  rheology:  for  example, 
the  frequency-independence  of  At  constrains  the  value  of  the  exponential  a  to  —  1  (while 
the  parameter  £  becomes  simply  At).  However,  for  an  arbitrary  a  ^  —1  the  lags  will  read 
(Efroimsky  and  Lainey  2007): 


At 


lmpq 


-a  — (a+1) 

=  £  x. 

lmpq 


(83) 


While  the  MacDonald  approach  cannot  be  generalised  to  a  ^  —  1 ,  the  Darwin-Kaula-Goldr- 
eich  method  can  be  well  combined  with  (83).  To  this  end,  we  shall  insert  (80-81)  and  (83) 
into  (79a),  and  shall  also  employ  the  evident  formulae 


|  cote 


cose  = 

lmpq 


lmpq 


Q 


lmpq 


f  x 


lmpq 


1  +  cot2  e 


lmpq 


1  +  Q- 


lmpq 

sgn  co 


t  .  c.2a  2a 

l+£  x, 

lmpq 


(84) 


sin  6  =  sin  16  1  sgn  co  = 

lmpq  lmpq  lmpq 


lmpq 


sgn  co 


lmpq 


1  +  cot2  e 


lmpq 


1  +  Q 


lmpq 


sgn  co 


lmpq 


1  .  ^-.2  a  2  a 

l+£  x, 

lmpq 


(85) 


with  tu  given  by  (25),  and  |e  ]  assumed  (for  reasons  explained  in  the  Appendix)  not  to 
approach  too  close  to  7r/2.  This  will  give  us  the  following  expression  for  (the  constant  part 
of)  the  deceleration  rate  of  a  non-resonant  prograde  spin: 
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£  \29  —  4n\ 


+  - e4sgn  (29  —  An)  ,  .  , 

4  6  1  +r“|20  -4«|-“ 

+  Od2/Q)+  0(e6/Q). 


(86) 


Be  mindful,  that  a  naive  substitution  of  the  formula  (85)  for  sin  elmpq  into  (79a)  would  result 
in  an  expression  for  the  torque,  attaining  its  maxima  on  approach  to  resonances  (for  a  positive 
a),  an  evidently  unphysical  behaviour.  As  explained  in  Sect.  9,  there  exists  a  profound  physi¬ 
cal  reason,  for  which  the  actual  multiplier  in  (79a)  must  be  not  sin  elmpq  but:  sin  elmpq  cos  elmpq . 
Mathematically,  the  presence  of  the  cosine  is  irrelevant  unless  x]mpq  and  2lmpq  approach  zero. 
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If,  however,  x  becomes  very  small  (i.e.,  if  we  approach  a  resonance),  it  is  this  long-omit¬ 
ted  (though  known  yet  to  Darwin  1879)  cosine  multiplier  that  saves  us  from  the  unphysical 
maxima — See  Sect.  9  above. 

Under  the  extra  assumptions22  of  |elmpq|  1  and  n  -C  6,  formula  (90)  simplifies  to 


*=c=/c|- 
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where  the  overall  factor  K.  is  given  by 


K.  = 


3n2  M~ .IctAl. 
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(87) 

(88) 


(89) 


an  expression  identical  to  (78),  except  that  it  contains  A t22QQ,  Q,200,  and  x2200  instead  of  At, 
Q,  and  x,  correspondingly. 

Were  a  equal  to  —1,  sum  (87)  would  coincide  with  (82),  provided  6  >  2n  (but  not 
otherwise!).  For  realistic  mantles  and  crusts,  though,  the  values  of  a  will,  as  pointed  above, 
reside  within  the  interval  0.2  —  0.4  (closer  to  0.2  for  partial  melts). 


11  Conclusions 


In  the  article,  thus  far,  we  have  provided  a  detailed  review  of  a  narrow  range  of  topics.  Our 
goal  was  to  punctiliously  spell  out  the  assumptions  that  often  remain  implicit,  and  to  bring 
to  light  those  steps  in  calculations,  which  are  often  omitted  as  “self-evident”. 

This  has  helped  us  to  demonstrate  that  MacDonald-style  formula  (57)  for  the  tidal  torque 
is  valid  only  in  the  zeroth  order  of  enj Qx,  while  its  time-average  is  valid  only  in  the  first 


22  The  smallness  of  |eta  |  enables  one  to  employ  (79b)  instead  of  (79a).  Then  the  multipliers 
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=  sgn  0)22Qti  A t22Q0  |^(1  +  a)x2200«x220</  J,  the  latter  approximation  being  legitimate  only  under 
the  condition  of  X220</  —  X2200  X2200 »  which  turns  out  to  be  equivalent  to  n  <$C  0.  For  example, 
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we  arrive  at  formula  (87). 
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order.  These  restrictions  mean  that  in  the  popular  expressions  for  tidal  despinning  rate  the 
terms  with  higher  powers  of  e  become  significant  only  for  large  eccentricities.  Their  sig¬ 
nificance  is  limited  even  further  by  the  error  0((9~3)  emerging  when  the  sine  of  the  phase 
lag  gets  approximated  by  the  inverse  quality  factor — see  formula  (59)  and  the  paragraph 
preceding  it. 

We  have  demonstrated  that  in  the  case,  when  the  inclinations  are  small  and  the  phase 
lags  of  the  tidal  harmonics  are  proportional  to  the  frequency,  the  Darwin-Kaula  expansion  is 
equivalent  to  a  corrected  version  of  the  MacDonald  formalism.  The  latter  method  rests  on  the 
assumption  of  existence  of  one  total  double  bulge.  The  necessary  correction  to  MacDonald’s 
approach  would  be  to  assert  (following  Singer  (1968))  that  the  phase  lag  of  this  integral  bulge 
is  not  constant,  but  is  proportional  to  the  instantaneous  synodal  frequency  2(v  —  9),  where 
v  and  9  are  the  true  anomaly  and  the  sidereal  angle.  Any  rheology  different  from  this  one 
will  violate  the  equivalence  of  the  Darwin-Kaula  and  MacDonald  descriptions.  It  remains 
unexplored  if  their  equivalence  is  violated  also  by  setting  the  inclination  high. 

We  have  demonstrated  that  no  “paradoxes”  ensue  from  the  frequency-dependence 
Q  ~  x“  ,  with  a  =  0.3  ±  0.1,  found  for  the  mantle. 

We  have  investigated  the  limitations  of  the  popular  formula  interconnecting  the  quality 
factor  Q  and  the  phase  lag  e.  It  turns  out  that  for  low  values  of  the  quality  factor  (much 
less  than  10),  the  customary  formula  Q  =  cot  |e|  should  be  substituted  with  a  much  more 
complicated  relation. 

Finally,  we  examined  two  derivations  of  the  popular  expressions  (72-74),  and  have  pointed 
out  that  these  expressions  have  limitations  related  to  the  frequency-dependence  of  the  quality 
factor.  First,  dependent  upon  the  values  of  e  and  Q,  the  high-order  terms  in  these  expressions 
may  become  significant  only  for  large  eccentricities.  Second,  the  expansion  of  the  deceler¬ 
ation  rate  in  even  powers  of  e  will  be  different  if  At  is  frequency-dependent  (which  is  the 
case  for  solid  materials).  These  two  circumstances  do  not  necessarily  disprove  any  major 
result  achieved  in  the  bodily-tide  theory.  However,  some  coefficients  may  now  have  to  be 
reconsidered. 

For  the  realistic  rheology  of  terrestrial  bodies,  the  despinning  rate,  in  the  absence  of  tidal 
locking,  is  given  by  our  formulae  (86-89). 
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Appendix.  The  lag  and  the  quality  factor:  is  the  formula  Q  =  cot  |e|  universal? 

The  interrelation  between  the  quality  factor  Q  and  the  phase  lag  e  is  long-known  to  be 

Q  =  cot  |e|,  (90) 

and  its  derivation  can  be  found  in  many  books.  In  Appendix  A2  of  Efroimsky  and  Lainey 
(2007),  that  derivation  is  reproduced,  with  several  details  that  are  normally  omitted  in  the 
literature.  Among  other  things,  we  pointed  out  that  the  interrelation  has  exactly  the  form  (90) 
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only  in  the  limit  of  small  lags.  For  large  phase  lags,  the  form  of  this  relation  will  change 
considerably. 

Since  in  Sect.  9  of  the  current  paper  we  address  the  case  of  large  lags,  it  would  be  worth 
reconsidering  the  derivation  presented  in  Efroimsky  and  Lainey  (2007),  and  correcting  a 
subtle  omission  made  there.  Before  writing  formulae,  let  us  recall  that,  at  each  frequency  x 
in  the  spectrum  of  the  deformation,  the  quality  factor  (divided  by  2  jr)  is  defined  as  the  peak 
energy  stored  in  the  system  divided  by  the  energy  damped  over  a  cycle  of  flexure: 


2  JZ  F peak  (X  ) 
Afi  cycle  (x) 


(91) 


where  A£cycie(x)  <  0  as  we  are  talking  about  energy  losses.23 

An  attempt  to  consider  large  lags  (all  the  way  up  to  |e|  =  jr/2 )  sets  the  values  of  Q/ln 
below  unity.  As  the  dissipated  energy  cannot  exceed  the  energy  stored  in  a  free  oscillator, 
the  question  becomes  whether  the  values  of  Q/ 2jt  can  be  that  small.  To  understand  that 
they  can,  recall  that  in  this  situation  we  are  considering  an  oscillator,  which  is  not  free  but 
is  driven  (and  is  overdamped).  The  quality  factor  being  much  less  than  unity  simply  implies 
that  the  eigenfrequencies  get  damped  away  during  less  than  one  oscillation.  Nonetheless, 
motion  goes  on  due  to  the  driving  force. 

Now  let  us  switch  to  the  specific  context  of  tides.  To  begin  with,  let  us  recall  that  the 
dissipation  rate  in  a  tidally  distorted  primary  is  well  approximated  by  the  work  that  the 
secondary  carries  out  to  deform  the  primary: 


E 


pV  ■  VW  d3x 


(92) 


p ,  V,  and  W  denoting  the  density,  velocity,  and  tidal  potential  inside  the  primary.  The 
expression  on  the  right-hand  side  can  be  transformed  by  means  of  the  formula 

p  V  ■  VW  =  V  •  (pVW)  -WV-Vp-WV-  (pV) 

=  V  •  {pVW)  -  WV  ■  Vp  +  W  — ,  (93) 

3 1 

where  the  W  V  -V p  and  dp/dt  terms  may  be  omitted  under  the  assumption  that  the  primary 
is  homogeneous  and  incompressible.  In  this  approximation,  the  attenuation  rate  becomes 
simply 


E 


V  •  ( pVW)d3x 


pWV  ■  n  dA, 


(94) 


ii  being  the  outward  normal  to  the  surface  of  the  primary,  and  dA  being  an  element  of 
the  surface  area.  It  is  now  clear  that,  under  the  said  assertions,  it  is  sufficient  to  take  into 
account  only  the  radial  elevation  rate,  not  the  horizontal  distortion.  This  way,  formula  (92), 
in  application  to  a  unit  mass,  will  get  simplified  to 


E  = 


(95) 


23  We  are  considering  flexure  in  the  linear  approximation.  Thus  at  each  frequency  x  the  appropriate  energy 
loss  over  a  cycle,  A  i^cycle  (X ) »  depends  solely  on  the  maximal  energy  stored  at  that  same  frequency,  fipeak  (X )  ■ 
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f  standing  for  the  vertical  displacement  (which  is,  of  course,  delayed  in  time,  compared  to 
W  ).  The  amount  of  energy  dissipated  over  a  time  interval  ( t0 ,  t)  is  then 


A  E 


(96) 


We  shall  consider  the  simple  case  of  an  equatorial  moon  on  a  circular  orbit.  At  each  point 
of  the  planet,  the  variable  part  of  the  tidal  potential  produced  by  this  moon  will  read 


W  =  W0  cos  xU  (97) 

the  tidal  frequency  being  given  by 

X  =  2\n-cop\,  (98) 

n  being  the  secondary’s  mean  motion,  and  cop  being  the  primary’s  spin  rate.  Let  g  denote  the 
surface  free-fall  acceleration.  An  element  of  the  planet’s  surface  lying  beneath  the  satellite’s 
trajectory  will  then  experience  a  vertical  elevation  of 

W 

?  =  h2  —  cos(xf  -  |e|),  (99) 

g 

I12  being  the  corresponding  Love  number24,  and  |e|  being  the  positive 25  phase  lag,  which 
for  the  principal  tidal  frequency  is  simply  the  double  geometric  angle  S  subtended  at  the 
primary’s  centre  between  the  directions  to  the  secondary  and  to  the  main  bulge: 

|e|  =  2  8.  (100) 


Accordingly,  the  vertical  velocity  of  this  element  of  the  planet’s  surface  will  amount  to 


W0  Wo 

u  =  t  =  ~  h2  X  —  sin(xf  -  |e|)  =  -  h2  x  —  (sin  /f  cos  |e|  -  cos  xt  sin  |e|). 
g  g 

(101) 


The  expression  for  the  velocity  has  such  a  simple  form  because  in  this  case  the  instantaneous 
frequency  /  is  constant.  The  satellite  generates  two  bulges — on  the  facing  and  opposite  sides 
of  the  planet — so  each  point  of  the  surface  is  uplifted  twice  through  a  cycle.  This  entails  the 
factor  of  two  in  the  expression  (98)  for  the  frequency.  The  phase  in  (100),  too,  is  doubled, 
though  the  necessity  of  this  is  less  evident, — see  Appendix  A1  to  Efroimsky  and  Lainey 
(2007). 


24  For  a  homogeneous  incompressible  body,  kn  =  (3/5)h2,  for  which  reason  (99)  and  the  subsequent 
equations  with  h2  can  equally  be  written  as  proportional  to  k2.  The  formulation  employing  k2  is  more 
fundamental,  as  it  can,  in  principle,  be  generalised  to  a  compressible  body  with  a  radially  changing  density. 
Indeed,  whatever  the  properties  of  the  primary  body  are,  the  dissipation  rate  in  it  is  equal  to  the  rate  of  change 
of  the  primary’s  spin  energy  plus  the  rate  of  change  of  the  orbital  energy  in  the  system.  Both  the  latter  and  the 
former  are  proportional  to  k2. 

25  Were  we  not  considering  the  simple  case  of  a  circular  orbit,  then,  rigorously  speaking,  the  expression  for 

W  would  read  not  as  W0  cos  /  l  but  as  W()  cosa>ti(la|f,  the  tidal  frequency  ct>ti(jaj  taking  both  positive  and 
negative  values,  and  the  physical  frequency  of  flexure  being  /  =  In^dal  I-  Accordingly.  Lire  expression  for  J 
would  contain  not  cos  (xt  —  |e|)  but  cos(o)ti(jajf  —  € ).  As  we  saw  in  Eq.  24,  the  sign  of  e  is  always  the  same 
as  that  of  tuti(ja| .  For  this  reason,  one  may  simply  deal  with  the  physical  frequency  /  =  I  and  with  the 

absolute  value  of  the  phase  lag,  |e|. 
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The  energy  dissipated  over  a  time  cycle  T  =  2tt/x  ,  per  unit  mass,  will,  in  neglect  of 
horizontal  displacements,  be 


A  E  .  = 

cycle 


g  J  dr 


rt=T 


cos /t  (sinx?  cos|e|  —  cos xt  sin|e|)r/t 


Jt= 0 


=  -  h  2  x  — 


3  W0  . 


1  rxt=ln 


g  dr 


sin  e  — 


X  J xt= o 


cos 1  xt  d(xt) 


,  W0  dW„ 

=  —  /Z2 - 7T  sin  \e\, 

g  dr 


(102) 


while  the  peak  energy  stored  in  the  system  during  the  cycle  will  read: 


f  = 
^peak 


[T/A  / 

3W^ 

( 

/  U  ( 

\dt  =  - 

-  -  h2  x  — 

J  |e|/x  V 

dr  ) 

V  g  / 

rt=T/4 

X  / 

cos  xt  (sin  xt 

cos  |e|  —  cos  xt 

Jt=\e\/x 

3  W0 

"cos  |e| 

rxt=n/l 

X  hi  - 

/  cos  xt 

g 

dr 

X 

Jxt=\e  1 

sin  \€\ 


rXt=n/2 


>Xt=\e  I 


cos2  xl  d(xt ) 


(103) 


In  the  appropriate  expression  in  Appendix  A 1  to  Efroimsky  and  Lainey  (2007),  the  lower  limit 
of  integration  was  erroneously  set  to  be  zero.  To  understand  that  in  reality  integration  over 
Xt  should  begin  from  |e|,  one  should  superimpose  the  plots  of  the  two  functions  involved, 
cosx?  and  sin(x?  —  |e|).  The  maximal  energy  gets  stored  in  the  system  after  integration 
through  the  entire  interval  over  which  both  functions  have  the  same  sign.  Hence  xf  =  I6 1 
as  the  lower  limit. 

Evaluation  of  the  integrals  entails: 


dW0 

"i 

1 

l71  , 

dr 

2  cosl6l 

_  2 

sin  |e| 

whence 


_A  E  1 

q-  1  _  “^cycle  _  1 


jt  sin|e| 


2tt  E, 


peak 


2n 


^cos|6|  -  2  (j  -  h)  sinlf 


tan  lei 


1  -  (f  -  |e|)  tan  |e| 


(104) 


(105) 


As  can  be  seen  from  (105),  both  the  product  sin  e  cos  e  and  the  appropriate  component  of 
the  torque  attain  their  maxima  when  Q  =  1  —  jr/4. 

Usually,  |e|  is  small,  and  we  arrive  at  the  customary  expression 

Q-1  =  tan|e|  +  0(e2).  (106) 
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In  the  opposite  situation,  when  Q  — >  0  and  |e|  ->  7t/2  ,  it  is  convenient  to  consider  the 
small  difference 

§  =  |-kl,  (107) 

in  terms  whereof  the  inverse  quality  factor  will  read: 

,  cot  £11  1 

QTX  = - - —  = - = - = - , 

1  £  cot  ^  tan£-f  z  —  arctan z  1.  z3  |\  _  3  z2  +  0(z4)j  ' 

(108) 

where  z  =  tan£  and,  accordingly,  £  =  arctan  z  =  z  —  j  z3  +  ^  z5  +  0(z 7).  The  latter 
may,  of  course,  be  rewritten  as 


or,  the  same,  as 


l  -  -  r  +  0(z4) 


=  3  Q 


z  =  o  g)1/3 


1  +  ~  z2  +  0(z4) 


(109) 


(110) 


While  the  zeroth  approximation  is  simply  z  =  (3 Q) 1//3  +  O(Q)  ,  the  first  iteration  gives: 


tan  |  =  z  =  (3  Q) 


1/3 


i  +  T(3e)2/3  +  o(e4/3) 


1  +  -qz  +  0(q«) 


(111) 


with  q  =  (3  2)1/3  playing  the  role  of  a  small  parameter. 

We  now  see  that  the  customary  relation  (106)  should  be  substituted,  for  large  lags,  i.e., 
for  small26  values  of  Q ,  with: 


tan  |e|  =  (3Q)-1/3 


1  -  7  (3  G)2/3  +  0(g4/3) 


(112) 


The  formulae  for  the  tidal  potential  and  torque  contain  a  multiplier  sin  e  cos  e ,  whose  absolute 
value  can,  for  our  purposes,  be  written  down  as 


sin  |e|  cos  |e|  =  cos  £  sin£  = 


tan 


|  Q  [l  +  j?2  +  CH<74)] 


l+tan2f  1  +  q2  [l  +  0(g2)] 


=  (30 


1/3 


1  -  7T  (3  Q)2/3  +  0(fi4/3) 


*]' 


whence 


sine  cose  =  ±  (3 Q) 


1/3 


i  -  -o  g)2/3  +  o(g4/3) 


(113) 


(114) 


an  expression  vanishing  for  Q  — >•  0.  Be  mindful  that  both  e22oo  and  the  appropriate  compo¬ 
nent  of  the  torque  change  their  sign  on  the  satellite  crossing  the  synchronous  orbit. 


26  The  afore-employed  expansion  of  arctan  z  is  valid  for  ]z|  <  1.  This  inequality,  along  with  (108),  entails: 
Q  =  z  —  arctan  z  <  1  —  tz/4. 
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